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Abstract. This work reviews the recent developments in the field of Generalized 
Parton Distributions (GPDs) and Deeply virtual Compton scattering in the valence 
region, which aim at extracting the quark structure of the nucleon. We discuss the 
constraints which the present generation of measurements provide on GPDs, and 
examine several state-of-the-art parameterizations of GPDs. Future directions in this 
active field are discussed. 
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1. Introduction: PDFs, FFs and GPDs 
1.1. PDFs, FFs and GPDs 

Electron scattering, or more generally lepton scattering, has always been a powerful 
tool to investigate the structure of the subatomic world. This is due to the structureless 
nature of leptons and the way these latter interact with the target, i.e. mainly through 
the electromagnetic interaction which is described by an extremely precise theory: 
Quantum Electro-Dynamics (QED). 

Since the late 60's and the advent of electron accelerators in the GeV range, nucleon 
structure has been investigated by two main classes of electron scattering processes: 
inclusive scattering eA^ — )■ e'X, also called "Deep Inelastic Scattering" (DIS), and elastic 
scattering eN — )■ e'A^', where e (e') stand for an electron, or more generally a lepton, 
{N') for a nucleon and X for an undefined final state. 

We recall that Hofstadter was awarded the Nobel Prize in 1961 "for his pioneering 
studies of electron scattering in atomic nuclei" which revealed that the proton explicitly 
appeared as an extended object and not as a pointlike particle. These measurements 
have shown that as the momentum transfer increases in the elastic scattering, the cross 
section sharply decreases compared to the electron scattering on a pointlike charge. 
In contrast, for the inelastic scattering process it was found firstly at SLAC that the 
cross section at large momentum transfers does not show the sharp fall-off as elastic 
scattering but shows a scaling behavior. This so-called Bjorken scaling put in evidence 
the presence of pointlike charged constituents within the nucleons, for which the 1990 
Nobel prize was awarded to Friedman, Kendall and Taylor. Through Feynman's parton 
model, these consituents were then identified with the "quarks" introduced earlier by 
Gell-Mann (who was awarded the Nobel Prize in 1969 for his "Eightfold way"), based 
on theoretical symmetry considerations. 

In a first approximation, electron scattering proceeds through a one-photon 
exchange (we will keep this approximation in the whole of this work) and is characterized 
by = —{pe — Pe)^ > 0, the squared four-momentum transfered to the nucleon by the 
electron. The virtuality of the spacelike virtual photon can be thought of as the 
resolution or the scale with which one probes the inner structure of the nucleon. 

At sufficiently high Q^, the quark structure of the nucleon can be "seen" and the 
DIS process can be depicted by Fig. [1] (left panel), where the incoming lepton interacts 
with a single quark of the nucleon via the exchange of a virtual photon. The signature 
of such a pointlike and elementary photon-quark process is the Q^-independence of 
the amplitude of the process, i.e. the absence of a scale in the process. DIS results 
accumulated for more than 40 years show that this picture, so-called "scaling", starts 
to be valid already at ^ 1 GeV^. 

The complex quark and gluon structure of the nucleon, governed by the theory of 
strong interactions. Quantum Chromo-Dynamics (QCD), in its non-perturbative regime 
is then absorbed in structure functions. This is the concept of QCD factorization where 
one separates a point-like, short-distance, "hard" subprocess, from the complex, long- 
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Figure 1. Left panel: Deep Inelastic Scattering: at high virtuality of the photon 
exchanged between the electron and the nucleon, the photon interacts with a single 
quark of the nucleon. The struck quark escapes the nucleon and hadronizes leaving the 
final hadronic state undetermined. Right panel: Illustration of the associated non-local 
diagonal matrix element (p|V'g(0)CV'g(y)b) accessed in DIS. 

distance, "soft" structure of the nucleon. The calculation of such soft matrix elements 
directly from the underlying theory amounts to solve QCD in its non-perturbative 
regime, which is still a daunting task. Ab initio calculations, by evaluating QCD 
numerically on a discretized space-time euclidean lattice, are at present the most 
promising avenue to provide predictions for some category of such non-perturbative 
objects. In the DIS process, the soft structure functions are the well-known unpolarized 
and polarized Parton Distribution Functions (PDFs) q{x) and Ag(a;) respectively. In a 
frame where the nucleon approaches the speed of light in a certain direction, x is the 
longitudinal momentum fraction carried by the quark which is struck by the virtual 
photon. The PDFs represent therefore the (longitudinal) momentum distribution of 
quarks in the nucleon. 

The PDF structure functions correspond to QCD operators depending on space- 
time coordinates. Precisely, the PDFs are obtained as one- dimensional Fourier 
transforms in the lightlike coordinate y~ (at zero values of the other coordinates) as : 



q{x) 
Ag(a;) 



ixp 



=y±=o 



{pS\\\i^q{oh+i5i^q{y)\pS\i) 



where ipg is the quark field of flavor q, p represents the initial (and final, since it is the 
same for DIS by virtue of the optical theorem) nucleon momentum, x is the momentum 
fraction of the struck quark and S\\ is the longitudinal nucleon spin projection. 

One uses here the light-front frame where the initial and final nucleons are coUinear 
along the 2;-axis and the light-cone components are defined by = (a° ±a^)/V2. Since 
the space-time coordinates of the initial and final quarks are different, the operator in 
Eq. ([1]) is non-local, and since the momenta of the initial and final nucleons are identical, 
it is diagonal. This operator is illustrated in Fig. [T] (right panel) . 
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Figure 2. Left panel: elastic scattering. At high virtuality of the photon exchanged 
between the electron and the nucleon, the photon interacts with a single quark which 
remains in the nucleon. The nucleon has changed its momentum in the process 
but it remains a nucleon, unlike the DIS process where it has been "smashed into 
pieces" . Right panel: illustration of the associated local non-diagonal matrix element 
{p'\4^q{Gi)Oiljq{Q)\p) accessed in elastic lepton-nucleon scattering. 



The elastic eA^ — )■ eA^ process is illustrated in Fig. [2] (left panel). For this process, 
the long-distance "soft" physics is factorized in the Form Factors (FF) F2it), 
G\{t) and Gp(t), where t = {p^- — v'nY — ■ Iii the light-front frame, the squared 
momentum transfer t is the conjugate variable of the impact parameter. In such frame, 
the FFs reflect, via a Fourier transform, the spatial distributions of quarks in the plane 
transverse to the nucleon direction [H 121 El H] , see Ref. [5] for a recent review. 

The FFs are related to the following vector and axial-vector QCD local operators 
in space-time coordinates: 

(p'|^,(0)7+^,(0)b) = Fi^(t)iV(p')7+iVb) + F|(t)iV(p')za+'^-^iV(p), 

(p'|^,(0)7+75^,(0)b) = G^(t)iV(p')7+75iV(p) 

+ G?,(t)iV(p')75^iV(p), (2) 

where N and iV are the initial and final nucleon spinors and ttiat the mass of the nucleon. 
Such operators are "local" since the initial and final quarks are created (or annihilated) 
at the same space-time point and "non-diagonal" since the momenta of the initial and 
final nucleons are different. These operators are illustrated on the right panel of Fig. |2j 

PDFs and FFs have been measured for the last 40 years but are still an intense 
subject of investigation. The behavior of PDFs at large x is still a mystery and the 
recent observation of a different Q^-dependence for the electric and magnetic FFs was 
a true surprise (see Refs. [SI E] for recent reviews). 

A new avenue in the study of nucleon structure has opened up over the 
past two decades with the investigation of exclusive electroproduction processes. 
Theoretically, the formalism of Generalized Parton Distributions has emerged in the 
90's and experimentally, the latest generation of high-intensity and high-energy lepton 
accelerators, combined with high resolution and large acceptance detectors, allows to 
access such exclusive processes in a precise and systematic way. 
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Deep Exclusive Scattering (DES), i.e. the exclusive electroproduction of a photon 
or meson on the nucleon at large Q^, is illustrated on the left panel of Fig. |3] for 
the case of Deeply Virtual Compton Scattering (DVCS). The theoretical formalism 
and the factorization theorems associated with these processes have been laid out in 
Refs. [H [9l [m [121 [131 [H] • The corresponding factorizing structure functions are the 
so-called Generalized Parton Distributions (GPDs) Hi{x,^,t), E«(x,^,t), ^«(x,^,t) 
and E'^{x,C,,t). They correspond to the Fourier transform of the QCD non-local and 
non-diagonal operators which are illustrated on the right panel of Fig. [3] : 

^ f dy-e''''"'y-{p\^lj,{Oh+ij,{y)\p) 

y+=y±=o 

= H''{x,U) iV(p')7+iV(p) + E'^{x,U) N{p)ia+''p^N{p) , 
^ [ dy-e^^''^y'{p\^,{Oh+^'ijM\p) 

^ y+=y^=0 

= H\x, e, t) N{p)^^^,N{p) + E\x, e, t) Nip'H^Nip) , (3) 

where P is the average nucleon 4- momentum: P = {p + p')/2 and A = p' — p, the 4- 
momentum transfer between the final and initial nucleons. The combination of variables 
X + ^ is the light-cone -|— momentum fraction (of P) carried by the initial quark and the 
combination x — ^ is the -|— momentum fraction carried by the final quark going back in 
the nucleon. The variable t, the squared 4-momentum transfer between the final nucleon 
and the initial one, is defined as A^. GPDs depend on additional variables compared 
to PDFs and FFs. They are therefore a richer source of nucleon structure information, 
which we will detail in the following subsection. 

The QCD operators of Eq. (jS]) are "non-local" since the initial and final quarks 
are created (or annihilated) at different same space-time points and "non-diagonal" 
since the momenta of the initial and final nucleons are different. These operators are 
illustrated on the right panel of Fig. |3l 




m 



N{p) V t = A= ^ 7V(p + A) 



Figure 3. Left panel: the "handbag" diagrams for DVCS. The factorization theorems 
state that this is the dominant process at sufficiently high virtuality of the initial 
virtual photon: it is the same quark that has been struck by the virtual photon and 
that radiates the final photon. Right panel: illustration of the associated non-local 
non-diagonal matrix element (p'|7/'^(0)C'?/'g(2/)|p) accessed in DES. 
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The leading DVCS amplitude in the hard scale Q, the so-called twist-2 amplitude, 
corresponds to the transition between transverse photons. The 7^ — )■ 7^ transition is 
of order 1/Q and involves higher-twist (twist-3) quantities. Such quantities will not be 
discussed in this review which focuses on a leading-twist description of DVCS. Most 
studies indeed rely on the twist-2 assumption, which allows a first interpretation of 
existing DVCS data as will be shown below. Moreover genuine twist-3 structures are 
rather poorly known, and the restricted range of present DVCS measurements does 
not allow a clean and simple separation of leading-twist and higher-twist contributions. 
Therefore twist-3 effects will not be discussed in this review although they are required 
to ensure QED gauge invariance [151 IlHl [13, [E] ■ More generally higher-twist effects are 
not taken into account. In particular we will not cover the recent results on target mass 
and finite-t corrections to DVCS [HI [20]. Even if these new results suggest potentially 
large corrections to the leading-twist DVCS amplitude, they have not be included in 
any phenomenological study yet and their discussion is beyond the scope of this paper. 

In this review, we also focus on the quark helicity conserving quantities, i.e. the 
operators between the quark spinors in Eqs. ([1]) and ([3]) corresponding with 7"^ or 7+7^ 
matrices. One generalization involves the use of the a^^ operator, allowing to define 
"transversity" PDFs and GPDs. We will also concentrate only on quark GPDs, as we 
are interested in the valence region in the present work. One can also define gluonic 
GPDs corresponding to the operators: 



are illustrated in Fig. [H They will not be considered further in this review which is 
devoted to the study of DVCS in the valence region. At leading-twist gluon GPDs 
contribute at next-to-leading order in the strong coupling constant as- It is commonly 
believed that they have a small impact in the valence region, hence justifying the leading- 
order approximation. However complete next-to-leading order calculations of DVCS 
are available [2ll[22l[23l[2ll[25l[26l[271[28l[29] and recent estimates [30] challenge 
the common view: gluon contributions may not be negligible even in the valence 
region at moderate energy. These results triggered an ongoing theoretical effort on 
the soft-coUinear resummation in DVCS [31] 132]. New developments in this direction 
are expected in the near future but it is too early to detail them further here. 

We summarize in Table [T]the quark operators and the associated structure functions 
that we have just discussed. 

We refer the reader to Refs. [33l [311 [SSI [33 [371 [3H] for complete reviews on the 
GPD formalism. 

1.2. Properties of GPDs 

In Eq. ([3]), the GPDs H and E correspond with averages over the quark helicity. They 
are therefore called unpolarized GPDs. The GPDs H and E involve differences of quark 
helicities and are called polarized GPDs. At the nucleon level, E and E are associated 



(p'|G+^'(0)G+(y)b) and (p'|G+^(0)G+(y)b) , 
where G^" is the gluon field tensor and G^'' = le'^'^^'^Ga/s 
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HJ,EJi{x,e,,t) 
H9,E3(x,l,t) 

t = A2 




N(p + A) 



Figure 4. Illustration 
(p'|G+^(0)G+(y)b). 



of the non-local non-forward matrix element 



Operator 
in coordinate space 


Nature of the 
matrix element 


Associated structure functions 
in momentum space 


{p\mOHy)\p) 


non-local, diagonal 


fi{x),gi{x) 


{p'\mom\p) 


local, non-diagonal 


Fi(t),F2(t),GA(t),Gp(t) 


{p'\moHy)\p) 


non-local, non-diagonal 


H{x, e, t), E{x, t), H{x, t), e{x, e, t) 



Table 1. The three families of operators discussed in this section (with O = 7+ or 
7V). 



to a flip of the nucleon spin while H and H leave it unchanged. The four GPDs 
therefore reflect the four independent helicity-spin combinations of the quark-nucleon 
system (conserving quark helicity). These are illustrated in Fig. |5l 

Omitting the Q^-dependence associated with QCD evolution equation, the GPDs 
depend on three independent variables: x, ^ and t. x varies between -1 and 1 and ^ 
in principle also between -1 and 1 but, due to time reversal invariance, the range of ^ 
is reduced between and 1. If \x\ > ^, GPDs represent the probability amplitude of 
finding a quark (or an antiquark if x < in the nucleon with a -|- momentum fraction 
X + ^ and of putting it back into the nucleon with a -|- momentum fraction x — ^ plus 
some transverse momentum "kick", which is represented by t (or A^^). 

The other region — ^ < x < ^ implies that one "leg" in Fig. [3] (right panel) has 
a positive momentum fraction (a quark) while the other one has a negative one (an 
antiquark). In this region, the GPDs behave like a meson distribution amplitude and 
can be interpreted as the probability amplitude of finding a quark- antiquark pair in 
the nucleon. This kind of information on qq configurations in the nucleon and, more 
generally, the correlations between quarks (or antiquarks) of different momenta, are 
relatively unknown, and reveal the richness and novelty of the GPDs. 

Each GPD is defined for a given quark flavor: H"^ , E"^ , H'^ , E'^ {q = u,d,s,...). H 
and H are a generalization of the PDFs of Eq. ([1]) measured in DIS. In the forward 
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H : ^ 






E : ^ 






: ^ 




^ — ^. 





E : ^ 





p p 





p' 



Figure 5. The four GPDs H , E, H and E correspond to the various quarks hehcity 
and nucleon spin orientations. 



direction, one has the model independent relations : 

q{x), X > 0, 

x), a; < . 



i/''(x,0,0) 
^''(x,0,0) 



(5) 
(6) 



Ag(x), a; > , 
Ag(— x), X < 0. 

The origin of these relations is the optical theorem and the symmetry of the forward 
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Compton process, corresponding with zero four-momentum transfer, i.e. ^ = and 
t = 0. Fig. [6] illustrates this relation. 



N 




OC Xm 




N 



N 



t=0 



Figure 6. The optical theorem : the cross section of the DIS process is equal to the 
imaginary part of the forward amplitude of the (doubly virtual) Compton process. 



Also, at finite momentum transfer, there are model independent sum rules which 
relate the first moments of GPDs to the elastic FFs of Eq. ([2]): 

dxH^ix, e, t) = Flit) , dxE\x, e, t) = Flit) , 



dxH^ix^^^t) = G\{t) 



Fig. [7] illustrates these sum rules. 



dxE''{x,^,t) = G'^(t) . 



(7) 




Figure 7. Illustration of the sum rules of Eq. ([7]) linking the first x moment of the 
GPDs to FFs. The integration on the x momentum leads to a local operator (via the 
S function in y). 

One therefore sees that the PDFs and the FFs appear as simple limits or moments 
of the GPDs. 

Similarly to the FFs, the t variable in GPDs is the conjugate variable of the impact 
parameter (in the light-front frame) [H |39l HQ]. For ^ = (where t = — A^), one has 
threfore an impact parameter version of GPDs through a Fourier integral in tranverse 
momentum Aj^ : 

H'^ix, b^) = I ^e-^-^-i/^(a:, 0, -A^^). (8) 



GPDs in the valence region from DVCS 



12 





DSSV08 [41J 


BBIO 


LSSIO [43J 


AAC08 


NNRR12 


15] 


AS(Q2) 


0.25 ±0.02 
-0.10 ±0.16 


0.19 ±0.08 
0.46 ± 0.43 


0.21 ±0.03 
0.32 ±0.19 


0.24 ±0.07 
0.63 ±0.81 


0.31 ±0.10 
-0.2 ± 1.4 



Table 2. Comparison of different determinations of the first moments of the singlet and 
gluon polarised PDFs at the scale = 4 GeV^ in the MS scheme. All uncertainties 
shown are statistical only. 



At .^=0, the GPD(x, 0, t) can then be interpreted as the probabihty of finding a parton 
with longitudinal momentum fraction x at a given transverse distance (relative to the 
transverse cm.) in the nucleon. In this way, the information contained in a traditional 
parton distribution, as measured in DIS, and the information contained within a form 
factor, as measured in elastic lepton-nucleon scattering, are combined and correlated in 
the GPD description. 

The second moment of the GPDs is relevant to the nucleon spin structure. It 
was shown in Ref.[9] that there exists a (color) gauge-invariant decomposition of the 
nucleon spin: 1 = Jg + Jg, where Jg and Jg are respectively the total quark and gluon 
contributions to the nucleon total angular momentum. The second moment of the 
GPD's gives (Ji's sum rule) : 

jg = ^ j^dxx [/f'?(x,e,t = o)±E«(x,e,t = o)]. (9) 

The total quark spin contribution Jg decomposes (in a gauge invariant way) as 
Jg = ^AS ± Lg where 1/2 AS and Lg are respectively the quark spin and quark 
orbital contributions to the nucleon spin. AS can be measured through polarized DIS 
experiments, and its extracted value is shown in Table HI One sees from Table H] that 
the different determinations of AS all point to a value in the range 20 - 30 %. 

On the other hand, for the gluons it is still an open question how to decompose 
the total angular momentum Jg into orbital angular momentum, Lg, and gluon spin, 
Ag, parts, in such a way that both can be related to observables. For a discussion on 
recent developments in this active field, see [16] and references therein. At present, it 
is only known how to directly access the gluon spin contribution Ag in experiment. It 
can be accessed in several ways: inclusive c — c or high pt hadron pairs production 
in polarized semi-inclusive DIS, semi-inclusive vr*^, 7, jet,... production in polarized 
proton collisions and evolution of gi{x,Q'^) through global fits of polarized data. At 
present Ag can only be extracted with a large uncertainty, as can be seen from Table [2j 
While most determinations for Ag indicate a too small value to fully explain the spin 
puzzle, it is clearly very worthwhile to reduce the uncertainty in its extraction by further 
measurements. 

The sum rule of Eq. ([9]) in terms of the GPDs provides a model independent way of 
determining the quark orbital contribution to the nucleon spin and therefore completes 
the quark sector of the "spin-puzzle" . 

Eq. (Q is actually a particular case of a more general rule on the x moments of 
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GPDs. The so-called polynomiality condition state that the x" moment of GPDs must 
be a polynomial in ^ of order n (for n even, corresponding with non-singlet GPDs) or 
n + 1 (for n odd, corresponding with singlet GPDs), e.g. for the H GPD : 

if n even : J dx x^H{x,$,,t) = Oq + a2^^ + a^^'^ + ... + an^"', 

if n odd: j dx x'^H{x,i,t) = + a2^'^ + a^i^ + ... + an+iC^^ . (10) 

There are similar rules for the GPDs i?, H and E. For the GPD i?, the a^+i coefficient 
is the same as for H except that it has the opposite sign. For the GPDs H and the 
maximum ^ power in Eq. ffTOl) for singlet GPDs is n — 1 (instead of n -|- 1). 

We note that in Eq. (ITOl) only even powers of ^ appear which is a consequence of 
the time reversal invariance which states that : H{x, —C,,t) = H{x,C,,t). 

Besides the polynomiality constraints on GPDs, the GPDs are also constrained 
by positivity conditions which should be taken into account both for nonzero and zero 
skewness parameter. The simplest of these conditions arises from requiring the positivity 
of the quark distribution in a transversely polarized nucleon. This imposes a relation 
between the ii^-type and if- type GPDs, for more details see Ref. [17] . 

As mentioned above, an ab initio calculation of soft matrix elements in general 
seems at present only pratical within lattice QCD. By its nature of discretizing the 
theory on an euclidean space-time lattice, lattice QCD can robustly calculate a few 
lowest moments of GPDs, which correspond to matrix elements of local operators. 
Although at present the calculations are being performed for unphysical pion masses, it 
is foreseeable that in the next few years calculations for such quantities at the physical 
point with controlled systematic uncertainties will become available. They can be used 
in the near future as additional constraints when confronting GPD parameterizations 
with experiment. We refer to the review paper of Ref. [H] for a recent review of lattice 
efforts in this field. 

1.3. Generalized Transverse- Momentum dependent parton Distributions (GTMDs) 

The GPDs can be considered as a particular limit of generalized parton correlation 
functions. The Generalized Parton Correlation Functions (GPCFs) provide a unified 
framework to describe the partonic information contained in a hadron. The GPCFs 
parameterize the fully unintegrated off-diagonal quark-quark correlator, depending on 
the full 4-momentum k of the quark and on the 4-momentum A which is transferred 
by the probe to the hadron; for a classification see refs. [l9l [50] . They have a direct 
connection with the Wigner distributions of the parton-hadron system [STl [521137] . which 
represent the quantum mechanical analogues of the classical phase-space distributions. 

When integrating the GPCFs over the light-cone energy component of the 
quark momentum one arrives at generalized transverse-momentum dependent parton 
distributions (GTMDs) which contain the most general one-body information of partons, 
corresponding to the full one-quark density matrix in momentum space. These GTMDs 
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parameterize the following general unintegrated, off-diagonal quark-quark correlator for 
a hadron [50] : 

12„ 



x{p',Al*(-|)rW«i(|)|p,A>|,,.„, (11) 

where the superscript V stands for any element of the basis {1, 75, 7^, 7^^755 
in Dirac space, and A (A') denote the helicities of initial (final) hadron respectively. 
A Wilson line W = VV(— |, ||n) ensures the color gauge invariance of the correlator, 
connecting the points — | and | via the intermediary points — | -|- 00 • n and | -|- 00 • n by 
straight lines. This induces a dependence of the Wilson line on the light-cone direction 
n. Furthermore, the parameter rj = sign(n°) gives the sign of the zeroth component of 
n, i.e. indicates whether the Wilson line is future-pointing [rj = +1) or past-pointing 
[rj = —1). Clearly, such correlators generalizes the GPD correlators introduced in 
Eq. (|3]) by allowing the quark operator to be also non-local in the transverse direction, 
i.e. besides the GPD arguments x, ^, and t, the GTMDs also depend on the quark 
transverse momentum k±. 

At leading twist, there are sixteen complex GTMDs, which are defined in terms of 
the independent polarization states of quarks and hadron. In the forward limit A = 0, 
they reduce to eight transverse-momentum dependent parton distributions (TMDs) 
which depend on the longitudinal momentum fraction x and transverse momentum 
kj_ of quarks, and therefore give access to the three-dimensional picture of the hadrons 
in momentum space. On the other hand, the integration over k± of the GTMDs leads to 
eight GPDs which are probability amplitudes related to the off-diagonal matrix elements 
of the parton density matrix in the longitudinal momentum space. After a Fourier 
transform of to the impact-parameter space, they provide a three-dimensional 
picture of the hadron in a mixed momentum-coordinate space. The common limit 
of TMDs and GPDs is given by the standard parton distribution functions (PDFs), 
related to the diagonal matrix elements of the longitudinal-momentum density matrix 
for different polarization states of quarks and hadron. Fig. [8] illustrates how the GTMDs 
reduce to different parton distributions and form factors. 

Although it has not been shown to date that the GTMDs can be accessed in a 
model independent way in experiment, they may however provide useful quantities to 
gain insight through model calculations of hadron structure, see e.g. [53] . 
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Figure 8. The GTMDs reduce to different parton distributions and form factors. By 
integrating over the quark transverse momentum fcj^, the GTMDs rciducc to the GPDs, 
whereas the forward hmit ^ = 0, f = (i.e. A = 0) can in turn be parameterized in 
terms of TMDs, which are the quantities which enter semi-inclusive deep inelastic 
scattering processes. Integrating the TMDs over k±_ gives rise to the forward parton 
distributions. The forward parton distributions arc also obtained by taking the forward 
limit A = in the GPDs. Integrating the GPDs over x yields the form factors. In the 
middle row right column figure, "fP stands for "fragmentation functions" which we do 
not discuss here. 
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2. Prom theory to data 

Among the hard exclusive leptoproduction processes, the DVCS channel bears the 
best promises to extract the GPDs from experimental data using the leading-twist 
handbag diagram amplitude. Indeed, in DVCS, the hard perturbative part of the 
handbag involves only electromagnetic vertices (Fig. [3]) while in Deeply Virtual Meson 
Electroproduction (DVMP), there are strong vertices involving a gluon exchange, see 
Fig. O In DVMP there is another soft non-perturbative quantity besides the GPDs that 
enters the calculation: the distribution amplitude (DA) of the meson which is produced. 
As the gluon virtuality needs to be hard to ensure the leading-twist amplitude, the 
endpoint behavior of the DA can potentially lead to strong power corrections to the 
leading-twist amplitude. In this review, we focus on the DVCS process on the proton. 




Figure 9. The leading-twist handbag diagram for DVMP involving quark (left) and 
gluon (right) GPDs. We note that the factorization between a hard scattering process 
and GPDs and DAs for DVMP has been demonstrated only for longitudinally polarized 
incoming photons. 



2.1. Compton Form Factors 

Four independent variables are needed to describe the 3-body final state reaction 
e{pe)p{Pp) ^' {p'e)pip'p)li.'P')) ^ fixed beam energy Ef,. They are usually chosen 
as Q^, xbi t and where = —{pe — Pe)^; = (with q the virtual photon four- 
momentum), t = {pp — p'yY and is the azimuthal angle between the electron scattering 
plane and the hadronic production plane (see Ref. [M] for its explicit definition within 
the Trento convention). See Fig. [10] for an illustration of these kinematical quantities. 

At leading- twist, the GPDs depend on the three variables : x, ^ and t where 
the variable ^ is related to xb by: ^ = xsji^ — xb)- In principle, GPDs depend on 
as well. However, this dependence can be predicted and calculated through the 
evolution equations and does not refiect the non-perturbative structure of the nucleon. 
For simplicity, we will not write the dependence on explicitly in the following. 
The variables ^ and t can be accessed by measuring the kinematics of the scattered 
electron and of the final state photon and/or proton. However, the variable x is not 
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experimentally accessible. In the calculation of the DVCS amplitude of the handbag 
diagram of Fig. [3l the variable x is integrated over. The DVCS amplitude is written as: 

Mdvcs = ^MyuiVi)Hl''o,Dvcs (12) 
where e^(g) and e*(p^) are respectively the polarisation 4- vectors of the (virtual) initial 
and final photons and: 



^L.O.DVCS 



2 J-i 



+ 



1 



X 



^ + ie x + ^ 



H 



^y^.ix, t) N{p h.riNip) + Ely^.ix, e, t) Nip )za^^^N{p) 



X — ^ + ie X + ^ — ie 



X 



^DVCS^ 



2m 



A ■ n 



X, e, t)N{p h.n^,N{p) + Elycsix, e, t)N{p )^,^N{p)i 1? 



2m 



with £0123 = +1 and the lightlike vectors along the positive and negative ^-directions: 
= P+/v^(l, 0, 0, 1) and = 1/P+ ■ 1/^2(1, 0, 0, -1). 

One readily sees from Eq. ( |T3l) that the variable x which is a "mute" variable is 
integrated over. It is also weighted by the ^^^-^^^ factors, which originate from the 
propagator of the quark in the handbag diagram of Fig. [3] (left panel). 

The DVCS amplitude contains convolution integrals of the form : 



1 x-^ + ie 



dx t7vH{^,^,t), 

1 



x-C. 



(14) 
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and analogously for the GPDs H or E. In Eq. fll4p . we have decomposed the 
expression into a real and an imaginary part where V denotes the principal value integral. 
This means that the maximum information that can be extracted from the experimental 
data in the DVCS process at a given t) point is GPD{±^, ^, t) or/ and J^^ dx ^^^^'^'^^ . 
The former is accessed when an observable sensitive to the imaginary part of the DVCS 
amplitude is measured, such as single beam- or target-spin observables, while the latter 
is accessed when an observable sensitive to the real part of the DVCS amplitude is 
measured, such as double beam- or target-spin observables or beam charge sensitive 
observables. The unpolarized cross section is sensitive to both the real and imaginary 
parts of the DVCS amplitude. 

There are therefore in principle eight GPD-related quantities that can be extracted 
from the DVCS process: 

HRei^, t)=V f dx [H{x, e, t) - H{~x, e, t)] C+{x, 0, (15) 
Jo 



Er,{U) =V [' dx[E{x,U) - E{-x,U)]C+{x,0, 
Jo 



(16) 







Hr,{U)=V dx H{x,U) + H{-x,U) C-{x,0, (17) 



ER,i^,t) = VJ dx[Eix,^,t) + Ei-x,^,t)^C-ix,0, (18) 
Hjmi^,t)^Hi^,^,t)-Hi-^,^,t), (19) 

EUU) =E{^,U)-E{-^,U), (20) 

HjU^, t) = iJ(e, t) + i^(-e, t), (21) 

EU^, t) = E(e, t) + Ei-^, t), (22) 
with the coefficient functions defined as : 

C^ix,0 = ^,±^ (23) 
X — t, X + t, 

and where one has reduced the x-range of integration from {—1, 1} to {0, 1} in the 
convolutions. 

The functions Hjie, Him, etc. on the Ihs of Eqs. (IT5]- [22|) . which depend on the two 
kinematical variables ^ and t, accessible in experiment, are called the Compton Form 
Factors (CFFs) Q . 

For further use, we also introduce the complex functions as : 

H(e,t) = HR,{^,t)-tnHjU^,t), (24) 

and analogously for the other GPDs. 

In the following, we will also use the notation: 

H+{x, e, t) = H{x, e, t) - H{-x, e, t), 

I Be aware of slightly different notations in the literature, e.g. the authors of Ref. [55] include — tt 
factors in the definition of the "/m" CFFs or include a minus sign in the definition of the "Re" CFFs. 
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E. 



E^ 



'+{x,^,t) =E{x,^,t)-E{-x,^,t), 
'+(x,e,t) =E{x,^,t) + E{~x,^,t), 



(25) 



(26) 



since these are the so-called singlet (C = +1) GPD combinations, which enter in the 
CFFs, and consequently in the DVCS observables. 

2.2. The Bethe-Heitler process 

The DVCS process is not the only amplitude contributing to the ep — )■ ejr^ reaction. 
There is also the Bethe-Heitler (BH) process, in which the final state photon is radiated 
by the incoming or scattered electron and not by the nucleon itself. This is illustrated 
in Fig. [121 The BH process leads to the same final state as the DVCS process and 
interferes with it. Since the nucleon form factors Fi and F2 can be considered as well- 
known at small t, the BH process is precisely calculable theoretically. The BH cross 
section has the very distinct feature to sharply rise around 0=0° and 180°. These are 
the regions where the radiated photon is emitted in the direction of the incoming electron 
or the scattered one. The strong enhancements when the outgoing photon is emitted in 
the electron-proton plane is due to singularities (for massless electrons) in the electron 
propagators. 

In those regions, a small variation in the kinematical variables Q^, xb-, t or 
produces strong variations in the cross section. For instance. Fig. [11] shows that a 
variation of xb by 1% around the particular kinematical setting xb=0.300, (^^=2.500 
GeV^, t=-0.200 GeV^, i?e=5.750 GeV induces a variation of almost 10% at very low 
and large 0. It should be noted that, even at 0=180°, this relative difference remains 
more than 5%. Experiments should therefore as much as possible quote and control the 
values of the kinematic variables at least at the percent level to achieve a few percent 
accuracy on the BH+DVCS cross sections. 

The importance of the BH relative to the DVCS strongly depends on the (x^, i-i 
0) phase space regions. It can largely dominate or be negligible compared to DVCS. The 
presence of BH can actually be considered as an asset when one measures observables 
which are sensitive to the BH-DVCS interference. The BH can then act as an amplifier 
for the DVCS process and gives access to CFFs in a linear fashion, instead of a bilinear 
one if only DVCS was present. 

2.3. Experimental observables 

As illustrated in Fig. [5l the four GPDs refiect the four independent spin/helicity 
nucleon/quark combinations in the handbag diagram. Therefore, the way to separate 
them is to use the spin degrees of freedom of the beam and of the target and measure 
various polarization observables. In Ref . [55] , analytical relations linking the observables 
and the CFFs have been derived for the ep — )■ ep7 reaction, considered as the (coherent) 



GPDs in the valence region from DVCS 



20 




Figure 11. The theoretical cross section of tlie BH process at the central kinematics 
Ee^b.lbQ GeV, .Tij^O.aOO, g2==2.500 GeV^, t=-0.200 GeV^, with variations of 1% in 
xb- 



DVCS 



BH 



a(ep —7- ep7) oc 




Figure 12. The cross section of the ep — >■ ep7 reaction is proportional to the squared 
amplitude : | Mjjycs + ^bh The DVCS process is proportional to GPDs while 
the BH process is proportional to FFs. The two processes interfere. 



sum of the BH and DVCS processes. We present here a few of them: 

Aaw oc sin0 Im{F{H + ^(Fi + F,)^ - kF2£ + ...} (27) 

AauL oc sin0 Im ^F{H + ^{F^ + F2) (^H + - ^A^Fa^ + ...} (28) 

AfTLL oc (A + 5COS0) i?e|Fi?i + e(Fi + F2) + + ...} (29) 

Aaux oc sin0 Im{k{F2'H - FiS) + ...} (30) 

where Act stands for a difference of polarized cross sections, with the first index 
referring to the polarization of the beam ( "U" for unpolarized and "L" for longitudinally 



GPDs in the valence region from DVCS 



21 



polarized) and the second one to the polarization of the target: "U" for unpolarized, 'L" 
for longitudinally polarized and "x" or "y" for a transversely polarized target. Indeed, in 
this latter case, there are two independent polarization directions: "x" is in the hadronic 
plane and "y" is perpendicular to it (see Fig. [T0|) . Furthermore, the kinematical variable 
k is defined as : k = — t/(4m^). 

The difference of polarized cross sections in Eqs. ( I271l30ll are sensitive only to 
the BH-DVCS interference term. The CFFs arise from the DVCS process while the 
Fi and F2 FFs originate from the BH process and only products of FFs and CFFs 
appear, thus providing access to CFFs in a linear fashion. At leading order in a ^ 
expansion, only sin or cos modulations appear. One also notices in general that single 
spin observables are sensitive to the imaginary CFFs while double-spin observables are 
sensitive to the real CFFs. 

In a first approximation, neglecting terms multiplied by kinematical factors such 
as ^, xb and fc, one can see that, on a proton target, I^oi^u is dominantly sensitive to 
E\^, AauL to ^f^, Actll to H^^^ and Adt/T to H^^^ and Ef^. For a neutron target, 
the sensitivity of these spin observables to the CFFs changes as the values of the FFs 
are different (in particular, Fi ^ at small t). Thus, on a neutron target, Actlc/ is 
dominantly sensitive to E]^^ and Hf^, ^ctul to Hf^^, ^cfll to H'^^ and /S.auT to Hf^. 

Here, we have displayed explicitely the "p" and "n" superscripts to underhne that 
GPDs on the proton and on the neutron are not equal. The relations expressing the 
"p" and "n" GPDs entering the DVCS amplitudes, in terms of the u- and (i-quark 
contributions, are given for the GPD H by : 



and similarly for the GPDs E, H or E. 

In summary, given the number of variables on which the GPDs depend (three, 
omitting the Q^-dependence), the convolution over x in the amplitudes, the presence 
of the BH process with it singularities, the number of CFFs (eight at leading twist), 
the quark flavor decomposition, not to mention the evolution and higher-twist 
corrections, it is clearly a non-trivial task to extract the GPDs from the experimental 
data and, ultimately to map them in the three variables x, ^, t. It requires a broad 
experimental program measuring several DVCS (or DVMP) spin observables on proton 
and neutron targets over large ranges in xb and t (and Q^). 

There are several strategies to make progress in such a program. One of them 
is, as an intermediate step, to extract the CFFs from DVCS data for a given (^, t) 
point by fitting the distribution at a given beam energy. This can be done 
in an essentially model-independent way provided one has enough constraints, i.e. 
experimental observables, to extract the eight CFFs. As we will see in section |U even 
if there are only two observables which are measured, to be fitted by eight CFFs taken 



H^i^, e, t) = ^H\i, e, t) + ]-H^{i, e, t), 



(31) 



(32) 
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as free parameters, making the problem a priori largely under-constrained, some results 
can still be obtained due to the domaince of such few observables by one or two CFFs. 
However, this is only the first step of the program, since the x dependence still needs 
to be uncovered, in principle with the help of a model with adjustable parameters. The 
problem can be simplified with the help of dispersion relations which we will discuss 
in section 13.41 They can in principle reduce from eight to five the number of GPD 
quantities to be extracted. They state, in a model-independent way, that the real 
part CFFs defined by Eqs. ([T5] - [18]) are actually integrals over ^ of their respective 
imaginary part CFFs defined by Eqs. (fTOl - [22l) . In this approach, there is in addition 
a real subtraction constant (at fixed ^ and t) which intervenes and makes the number 
of independent quantities to be five in total. To apply dispersion relations, it is needed 
to measure data over a very wide range in (at fixed t) unless one has good reasons 
to truncate the integral or to extrapolate. Another strategy consists in fitting directly 
the experimental observables by a model which has for each GPD H, E, H or E, a. 
parameterization of the full x, ^, t-dependence with parameters to be fitted. We will 
discuss these various approaches below. 

Let us also mention that there is an experimental way to measure indepedently the 
X and ^-dependence of GPDs. The double-DVCS process consists of the DVCS process 
with a virtual (space-like or time-like) photon in the final state. In the case of a final 
timelike photon, the virtuality of this second photon can be measured and varied, thus 
providing an extra lever arm and allowing to measure the GPDs for each x, ^, t values 
independently (though with some limitations if the final photon is timelike) [56l [57] . 
However, since the cross section of such process is reduced by a factor a ~ 1/137, and 
since one needs to make measurements above the vector meson resonance region to avoid 
the strong vector meson processes, the double DVCS has not revealed so far to be a 
practical way to access GPDs. 

We now review the existing DVCS measurements, limiting ourselves to the large 
and intermediate xb regions. 

2.4- Existing DVCS measurements 

Three experiments have provided these past 10 years DVCS data which can potentially 
lend themselves to a GPD interpretation. These are the Hall A and CLAS experiments 
from JLab (with a ~ 6 GeV electron beam energy) and the HERMES experiment at 
DESY (with a ~ 27 GeV electron or positron beam energy). 

2.4-1 ■ JLab Hall A The ep — )■ eV7 reaction was measured in the JLab Hall A 
experiment [58] by detecting only the scattered electron in a high resolution (y ~ 10"^ 
for momentum) arm spectrometer and the real photon in an electromagnetic calorimeter 
(.^ ^ 4 % for energy). A cut on the missing mass of the proton which clearly stood 
out over a small background was used to unambiguously identify the exclusive process. 
The Hall A experiment measured the 4-fold beam-polarized and unpolarized 
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differential cross sections da / dxBdQ'^dtdcf), i.e. without any integration over an 
independent variable, as a function of 0, for four —t values (0.17, 0.23, 0.28 and 0.33) 
at the average kinematics: < xb >= 0.36 and < >= 2.3 GeV^. The beam-polarized 
cross sections have also been measured at < >= 1.5 GeV^ and < >= 1.9 GeV^. 
Fig. [13] shows these results. The particular shape in of the BH contribution in the 
unpolarized cross section (red curve in the upper panels of Fig. [T3ll is easily recognizable. 
The difference between the red curve and the data is the contribution of the DVCS 
process and therefore of the GPDs. 



E^=5.75 GeV, Xg=0.36 




Figure 13. The panels of the top row show the DVCS/BH unpolarized cross 
sections on the proton, as a function of the (j) angle, measured by the JLab Hall 
A collaboration [58]. The average kinematics is < xb >=0.36, < >=2.3 GeV^ 
and < -t >=0.33, 0.28, 0.23 and 0.17 GeV^ (left to right). The red curves show 
the BH contribution. The bottom panels show the difference of beam-polarized cross 
section as a function of 4> for the same kinematics at < >=2.3 GeV^ (second row), 
< Q2 >^1.9 GeV^ (third row) and < >^i.5 Q^y^ (fourth row). 
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Figure 14. Top panel: Q^-dependence of the combination of CFFs: FiHj,n + £,{Fi + 
F2)Him ~ kF2Ejm for four different values of —t, as extracted from the fit of the beam- 
polarized cross sections of Fig. [13] (three bottom panels) . Bottom panel: i-dependence 
of the same combination of CFFs as extracted from the fit of the beam-polarized cross 
sections of Fig. [13] (three bottom panels). 



The difference of beam-polarized cross sections, i.e. Aaiu, is displayed in the three 
lower panels of Fig. [131 In the leading order ^ expansion, the amplitude of the sinusoidal 
is directly proportional to the combination of CFFs: FiT-L + C,{Fi + F2)'H — kF2£ (see 
Eq. (1271) ). Fitting these sinusoids has thus permitted to extract the Q^-dependence of 
this combination of CFFs at four different t values. The results of these fits is presented 
in Fig. [H] At leading- twist, GPDs and CFFs are predicted to be independent and 
the data seem to exhibit this scaling feature. Although the lever arm is very limited 
(~ 1 GeV^), this is a very encouraging sign that one can access the leading twist handbag 
process at the JLab kinematics. 

We also mention that the beam spin asymmetry of the DVCS-I-BH process on 
the neutron en — )■ en7 has been measured in an exploratory way by the JLab Hall A 
collaboration at one single (< xb >,< >) value (0.36,1.9) as a function of t [59] . 
Although these results are encouraging and might possibly give some first constraints 
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on the Eirn CFF, we decide, in this review, to focus on the proton channel. There is a 
variety of observables which have been measured over a wide phase space for this latter 
process. This should give the strongest constraints on the GPD models and fits. 

2.4-2. JLab Hall B The JLab CLAS collaboration uses a large acceptance spectrometer 
and has measured the DVCS process by detecting the three particles of the final state, 
i.e. the scattered electron, the recoil proton and the produced real photon, over a much 
broader phase space than in Hall A. Since CLAS has a lesser resolution ^ 10^^ 
for momentum) than the Hall A arm spectrometers, the kinematic redundancy and 
overconstraint due to the detection of the full final state is the best way to ensure the 
exclusivity of the process. 

Beam-polarized and unpolarized cross-sections measurements are under way [60] 
but to this day, only beam spin asymmetries, i.e. the ratio of Aam to the unpolarized 
cross section, and longitudinally polarized target asymmeties, i.e. the ratio of Aaui to 
the unpolarized cross section, have been measured. These asymmetries are observables 
which are relatively straightforward to extract experimentally since, in a first order 
approximation, normalization factors such as the efficiency/acceptance of the detector 
and, more generally, many sources of systematic errors cancel in the ratio. Both 
asymmetries have a shape close to a sin0 like Eq. ( l27|) and ( l28|) were predicting. The 
beam spin asymmetry was fitted by a function of the form a sin 0/(1 + c cos (p + d cos 20) . 
Fig. [15] (left panel) shows the value of this fitted asymmetry at = 90° for the ;^60 
{xBiQ'^it) bins covered by CLAS and which were measured with a 5.77 GeV beam. 

The longitudinally polarized target asymmetries are displayed in Fig. [161 These 
observables show a sin0-like shape, as predicted by theory (see Eq. (I28l)). and their 
sin moment A^^j^ is presented in the figure. Here we extend the subscript notation of 
Eqs. ( 1271) to (!30l) to asymmetry moments A, where the upperscript denotes the particular 
azumuthal moment considered. The use of a polarized target limited the statistics and 
the moments could be extracted only for 3 {xB,Q^,t) bins. 

2.4.3. HERMES At higher energies, xb ^ 0.1, the HERMES collaboration has carried 
out a measurement of ALL independent DVCS observables, except for cross sections: 
beam spin asymmetries [M| [65] , longitudinally polarized target asymmetries [66] . 
transversally polarized target asymmetries [671 [68], beam charge asjTiimetries [69l[70l[7T] 
and all associated beam spin/target spin and spin/beam-charge double asymmetries. In 
a first stage, the HERMES experiment requested the detection of the scattered electron 
(or positron) and of the final real photon. Then, a cut on the missing mass of the 
proton was applied. The width of this missing mass peak being more than 1 GeV, 
a substantial background subtraction had to be performed. In a second stage, the 
HERMES spectrometer was completed by a recoil detector allowing the detection of 
the recoil proton and therefore a complete identification of the DVCS final state. The 
kinematics being then over constrained, this allowed for a much cleaner selection of the 
exclusive reaction with a reduction of the contamination of non-DVCS events at the 
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Figure 15. DVCS-BH beam spin asymmetry at (f) — 90° as a function of t for different 
{xb,Q^) bins, as measured by the JLab Hall B/CLAS collaboration 61 (black solid 
circles) . The red empty triangles are the beam spin asymmetries derived from the ratio 
of the beam-polarized and unpolarized cross sections of Hall A (see Fig. [T5|) . The blue 
square point is the pionner measurement from the CLAS collaboration [52] , 

level of less than 1% [72]. In this "pure" DVCS samples, the amplitudes of beam-spin 
asymmetries were actually found to be somewhat larger (by about 10% in average). 

HERMES used a positron beam as well as an electron beam and the target spin 
asymmetries have a different sensitivity to the BH+DVCS amplitude according to the 
charge of the beam. To describe such correlated charge and beam-spin asymmetries, one 
therefore introduces a second index ( "I" or "DVCS" ) for the labelling of the asymmetries, 
according to (for instance for the beam-spin asymmetries): 
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Figure 16. The sin((/)) moments of the azimuthal moment of the DVCS-BH 
longitudinahy polarized target asymmetry [63]. Left: Three bins —t, integrated over 
xb] Right: Three bins in xb, integrated over t. 



( 4(0)-a^(0)) + (c^T(0)-^-(0) ) 
a+(0) + al(0) + a;(0) + orZ(0) 



^{LU,DVCS} — — , _+f±\ I _-f±\ I _-f±\ — ' I'^'^J 



^'""'^ a|(0) + al(0) + a;(0) + aZ(0) ' ^'^^ 

where the superscript represents the charge of the beam and the subscript the beam (or 
target) spin projection. At leading-twist, only the asymmetries with an "I" subscript 
can be sensitive to GPDs while the ones with the "DVCS" subscript are null. 

All DVCS azimuthal asymmetries have at leading twist a general sine, cosine or 
constant shape (see Eqs. [2711301 for a few examples), slightly modulated by the 0- 
dependence of the denominator. The sine, cosine and constant moments of the nine 
asymmetries which are expected to be non-null in the leading-twist handbag formalism 
are displayed in Fig. [T71 We added as a tenth observable the ^^i'^ moment (bottom 
right plot) which is expected to be power suppressed in this approximation. However 
the data show a 2- to 3-standard deviation difference from zero. If one doesn't consider 
this difference as a statistical fluctuation, it is a puzzle as it cannot be described by 
any leading-twist handbag calculation. Indeed, HERMES extracted also the "DVCS- 
subscript" asymmetries (see Eg. [Mil , as well as several sin 20, cos 20 or cos 30 moments, 
which are expected to be null in the hypothesis of DVCS leading-twist dominance. We 
do not display here these data but they were all found to be compatible with zero within 
error bars. Except for this puzzling A^^^^'^ moment, this gives further support to the 
idea that higher-twist contributions are small at the currently finite values explored. 
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Figure 17. Selection of ten DVCS-BH asymmetry i/i-moments as a function of t as 
measured by HERMES. All moments in this figure are expected to be non-null in the 
leading ^ expansion but for the A^^j^"^ moment (bottom right plot). 



confirming the first conclusions drawn from tlie JLab Hall A data. 

In Fig. [T71 we display only the t-dependence of these moments at the average x b 
and values of 0.09 and 2.5 GeV^ respectively. The data were taken with a 27.6 GeV 
beam energy. HERMES also measured the xb- and Q^-dependences, with the other 
kinematic variables fixed. Also, in this figure, the data correspond to data analysis 
carried out without the recoil detector. 

We also mention that the HERMES collaboration has measured the DVCS+BH 
charge, beam spin and longitudinally polarized target asymmetries with a deuterium 
target [731 E]- Such process is dominated by the incoherent DVCS+BH process on 
the proton and the results are in general consistent (with larger uncertainties) with the 
proton data shown in Fig. [T71 

We finish this section by mentioning that the unpolarized ep — )■ ejrf cross section 
has also been measured at much higher energy (30 < W < 120 GeV, 2 < < GeV^ 
where W is the center of mass energy of the 7* — p system), by the HI and ZEUS 
collaborations [751 US]- At such large W {i.e. low xb), the DVCS process is sensitive 
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mostly to "gluon" GPDs which we do not cover in this review. 
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3. Models of GPDs and dispersive framework for DVCS 

In this section, we review a few current state-of-tlie-art parameterizations of GPDs. We 
distinguish three famihes of models: models based on double distributions (VGG and 
GK), the dual parameterization and the Mellin-Barnes model. 

3.1. Double distributions / Regge phenomenology: the VGG and the GK models 

3.1.1. (x,C,) dependence and Double Distributions Double Distributions (DDs) were 
originally introduced by A. Radyushkin [771[78] and D. MuUer et al. [8]. They provide an 
elegant guideline to parameterize the [x,^) dependence of the GPDs which automatically 
satisfies the polynomiality relations (see Eq. f lTOj) ). 

The idea of the DDs is to decorrelate the transferred longitudinal momentum (A) 
from the initial nucleon momentum P (see Fig. [3]-right). In the light cone frame, one 
introduces then the new variables a and (3 such that the initial quark has a longitudinal 
momentum /3P+ - |(1 + a)A+ (see Fig. [IBl-left), instead of (x + ^)P+ (see Fig. |3]-left). 
Since, by definition, —2^ = this means that x = /3 + a^. The variable a is playing 
the role of ^, the fraction of longitudinal momentum of the transfer, but the difference is 
that a is now an absolute value, i.e. it has no reference to the (average) initial nucleon 
momentum, unlike ^. The link between a GPD and a DD is then only a change of 
variables, i.e. from (a, /3) to (x,^), such that : 

GPD\x,Cl= df3 da6{x- /3-^a)DD{a,f3) (35) 

J-i J-i+\P\ 



a 
1 




-1 



Figure 18. Left: illustration of the {p' \ {i;q{0)Otpq{y)\p) matrix element in terms of the 
a and /? variables, instead of x and ^ (see Fig. [3]-right). Right: the integration over 
the (a, /?) variables, in order to relate DDs and GPDs, takes place over the (3 = x — £,a 
straight line "inside" the rhombus defined by the equation | /3 | + | a | < 1- 
The slope of the straight line is given by ^ and its intercept at a = is given by x. 

One should integrate on all values/combinations of a and (3 which produce the (x, ^) 
variables. Given that x = f3 + ^a, one has actually only a one-dimensionnal integral. 
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The limits of the integration on the a and /3 variables are constrained by the fact that 
X has to be comprised between -1 and 1 and ^ between and 1, so that one has always 
I 3; I + I a |< 1. This constraint means that the integration over the variables a and 
(3 takes place over the {3 = x — straight line "inside" the rhombus defined by the 
equation | /? | + | a |< 1 (see Fig. [TSl-right). 

Due to the linear relation between x and ^ imposed by the S function, the 
polynomiality relation is automatically satisfied: the moment of Eq. ( 135|) will always 
produce a power. 

An advantage of the DDs is that the (a,/3) dependence can be more conveniently 
infered than the (x,^)-dependence. The matrix element corresponding to Fig. [TSl-left 
can be written: 

{p + A\^,{0)OMy)\p) (36) 

where we can consider two "extreme" cases. When there is no longitudinal momentum 
transfer brought by the photon, i.e. A = 0, one gets : 



{p\^g{OWg{y)\p) 



(37) 

y+=y±=o 



and one recovers, as it should, the forward matrix element of Eq. ([T]) and which is equal 
to the standard inclusive PDF, the forward limit of GPDs. 

However, there is now a second, new, limiting case : when P = and A 7^ 0, which 
is a case that couldn't be considered before, since A was proportional to P : 



{A\^P,{0)OMy)\0) 



(3^ 

y+=y±=o 



This matrix element should be interpreted as the probability amplitude to find in the 
nucleon a qq pair which shares the momentum A in (1 + a) and (1 — a) fractions (see 
Fig. [m with P = 0). Then, the idea is that the a functionnal dependence of the 
GPD could, in this domain, take the shape of a Distribution Amplitude (DA). A DA 
gives the probability amplitude to find in a meson M a qq pair which carries z and 
1 — z = z fractions of the meson momentum pm- The corresponding matrix element is 
the following : 



{0\mOij{y)\pM) 



(39) 

y+=y±=o 



where O can be a vector (7^) or axial (7'^7^) operator according to the parity of the 
meson. Such matrix element is illustrated In Fig. [19] and its Fourier transform reads : 

$m(^) = [ dy~e'^^P^'^y~ {0^0) Oij{y)\pM) • (40) 

y+=y±=o 

A DD can therefore be considered as a "mixture" /"hybrid" of a PDF and a DA, i.e. 
two limiting cases of a DD (respectively, A = and P = 0). Knowing the two limiting 
cases of the DD, the idea is to find a functional form of a and f3 which smoothly 
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Yl — <^ P 



Figure 19. Illustration of the (O|'0(O) O 



matrix clement. A DA 



-=y±=o 

represents the probability amplitude to find a quark and antiquark of momentum 
fraction, respectively, z and z — 1 — z, in a meson M of momentumpAf (or, equivalently, 
to create from the vacuum (0| a meson M with such a quark- antiquark pair). 



interpolates between a DA and a PDF when, respectively, a ^ and /9 — ?■ 0. A form 
which fulfills these requirements and proposed by Radyushkin [771 EH] , is : 

DD(/3,«) = /i(/3,a)g(/3), (41) 

h(R ^ r(26 + 2) [{l-\P\f-a^f 

'^^P^^) 22^+ir2(26 + l) (1-|/3|)2W ' 
where 6 is a free parameter. It governs the amount of ^ dependence of the DDs. The 
higher the b value, the weaker the ^ dependence for GPD'^{x,^,t). For instance, when 
b — )■ oo, h{P, a) — ?■ 1 and the DDs are independent of ^ and resemble a PDF. In principle, 
one can define a value for the valence, byai, and another one for the sea, bg^a- Figure [2U1 
shows H^{x, ^, t) as a function of x and ^ for t = and for b^ai = bgea = 1, following the 
DD ansatz of Eq. (135|) and (142|) . based on the VGG model which will be soon discussed. 

One identifies at ^ = a standard quark density distribution, with the rise around 
X = corresponding to the diverging sea contribution. The negative x part is related to 
antiquarks. One sees the evolution with ^ which tend towards the shape of an asymptotic 
DA. 

3.1.2. The D-term As we saw, GPDs built on the DD ansatz automatically satisfy 
the polynomiality rule. However, because of the 6{x — (3 — ^a) function in Eq. fl35|) . the 
^th ^ moment of the so defined GPDs is at most a polynomial in ^ of order n, while 
the polynomiality rule allows for a term with one more power, i.e. a ^"'''^ term. This 
means that the DD decomposition of the GPDs is not complete. The so-called D-teim, 
denoted by D{x/^, t) has been introduced by C. Weiss and M. Polyakov [79] to take into 
account this "missing" power It can be decomposed in a Gegenbauer series as : 

oo 

D{z,t) = {l-z')Y, dn{t)Cl/\z), (43) 

ri = l 

n odd 

with \z\ < 1. Since the D-term corresponds to a fiavor singlet contribution, it 
receives the same contribution from each quark fiavor. One can then define a D-term 
contribution for each quark fiavor by dividing by a factor Nj = 3, denoting the number 
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X 



Figure 20. The GPD H^{x, ^, t) as a function of the longitudinal momentum fraction 
X and the longitudinal momentum transfer ^ at i = according to the VGG model. One 
recognizes for ^=0 the typical shape of a parton distribution (with the sea quarks rising 
as X goes to 0, the negative x part being interpreted as the antiquark contribution) 
and as ^ increases the (asymptotic) shape of a distribution amplitude. 

of light quark flavors. Furthermore in Eq. (143|) . the form factors di{t), 6/3 (if:), ... at t = 
have been estimated, in a first approach, in the chiral soliton model as [34j : di = —4, 
c?3 = -1.2 et 4 = -0.4. 

The D-term "lives" only in the — ^ < x < ^ region, i.e. the qq part of the GPDs, 
whence the motivation to expand it on odd Gegenbauer polynomials which are the 
standard functions on which meson DAs are decomposed. We will come back to the 
D-term in section 13.41 devoted to dispersion relations. 

We are now going to describe the VGG and GK models which are both based on 
the DD (+ D-term) ansatz for the (x, ^)-dependence and which differ essentially by the 
parameterization of their t-dependence. 

3.1.3. The VGG model The VGG model is associated with a series of publications 
released between 1999 and 2005 [801 ED IMl El]. The first version of the model 
was published in 1998 but has since then continuously evolved, benefitting from and 
integrating the work and inputs of several other authors, as the field of GPD grew and 
improved. We quote here the main publications and associated steps and improvements 
of the VGG model over the past 10 years : 

• In Ref. [80], at first a parameterization of the GPDs via a ^-independent and 
t-factorized ansatz was used. In a concise notation : if (x, ^,t) = g(x)-Ff(t), 
Hi{x,^,t) = Aq{x)G\{t), etc... 

• In Ref. [81], the ^-dependence via the Double Distributions (DD) which was 
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discussed in the previous subsection, was introduced in the GPD parameterization. 
In a concise notation : Hi{x,^,t) = H%jj{x,OF^{t), Hi{x,^,t) = Hl^{x,OGA{t), 
etc... 

• In Ref. [34], the "D-term" (D), as suggested by C. Weiss and M. Polyakov [79], 
was introduced as well as a "Regge-inspired" unfactorized t-dependence, which are 
we are going to detail in the remaining of this subsection. In a concise notation : 
H'i{x,U) = Hl^{x,Ox-'''^' + ^Dil), E%x,U) = Elj,{x,i)x-'^' - 

etc... 

• In Ref. ^2], the Regge dependence was modified so as to satisfy the FF counting 
rules at large t. In a concise notation : H'^{x, ^, t) = ,^)a;~"'i^^~^'^* + 
Ei{x,^,t) = El^{x,0{^ - a;)^x-"^(i-)* _ etc... 

One can also complement this list by the work of Ref. [17] in which twist-3 effects 
are estimated in the Wandura-Wilczek approximation. 

(i) Parameterization of the GPD H 

The t-dependence of the GPD H in the VGG model is based on Regge theory. In 
short, Regge theory is based on the general concepts of unitarity and analyticity 
of scattering amplitudes and states that the high energy behavior of amplitudes 
should follow a s"''*^ behavior, where s is the (squared) center-of-mass energy of the 
system and a{t) is a Regge trajectory. A Regge trajectory is the relation between 
the squared mass and the spin of a family of mesons (or baryons) which share the 
same quantum numbers, except for spin. This high energy property is used, for 
instance, to determine, by extrapolation, the PDFs into the very small x region 
(i.e. high s domain since x ~ ■^) where no measurement is possible. 
Thus, q{x), which governs the x behavior of the DIS cross section, or, equivalently, 
of the forward {t = 0) Compton amphtude, should follow, as x — )■ 0, a -^^^ 
behavior. The leading meson trajectory associated to the valence part qv{x), which 
corresponds to an isovector combination {i.e. non-singlet), is the isovector vector 
meson p trajectory whose intercept ap{0) = 0.5. The sea (and gluon) part of the 
PDF qs{x) is an isoscalar (singlet) combination, and the associated trajectory is 
the Pomeron trajectory, which has the quantum numbers of the vacuum and for 
which ttp(O) = 1. One therefore infers that qv{x) ~ a;~°'^ and qs{x) ^ g{x) ^ x~^ 
as x — > (the trajectories governing the small x behavior of the polarized quark 
distributions are those of the axial vector mesons). 

For GPDs, i. e. the "non-forward" PDFs, the idea is to generalize this Regge ansatz 
for non-zero t values. Thus, the first formula which naturally suggests itself, is (for 
the H GPD and for ^ = 0, to simplify matter in a first stage) : 

H''{x,0,t) =q^{x)x-"'' (44) 

with the assumption of a linear Regge trajectory, i.e. : a(t) = a{0) + a't. a' is a 
parameter which can be strongly constrained by the sum rules linking the GPDs 
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to the Fi (and F2) FF(s), following Eq. ([7]) and, in particular, the nucleon charge 
radius. 

However, the ansatz of Eq. (jHj) has the shortcoming that it doesn't produce a 
correct behavior of the FFs at large t. At large t, quark counting rules dictate that 
Flit) should behave as ^. F2(t), which is spin flip, and is therefore suppressed, 
should behave as ^. The ansatz of Eq. ( Hll) does not satisfy these limits. Indeed, 
if Qvix) ~ (1 — xY for a; — > 1, one can show that, at large \t\: 



/•I 1 

/ (1 - xYx^'^'^dx oc — \t 
Jo « 



(45) 



which, with z/ ^ 3, taken from phenomenology, yields a asymptotic behavior 
for Fi{t), which is at variance with the behavior seen in the Fi form factor 
data, and expected from the asymptotic behavior. 

The large-t power behavior of Fi{t) should be governed by the large x (— )■ 1) 
behavior of q{x). Physically, the asymptotic large t-domain consists in probing the 
simplest configuration of the nucleon, i.e. the 3-quark configuration of the nucleon 
wave function which corresponds to the large x region. The idea is then to modify 
the large x behavior of Eq. (jH]). It can be done by introducing a (1 — x) term in 
the exponent of Eq. (jSj). With: 

\l - xYx-'^'^'-'-'^'dx oc - I t (46) 

a' 

one has, for z/ ^ 3, a ^ behavior at large t. 

To summarize, the ansatz for the VGG GPD H is therefore (for ^ = 0) : 

H{x, 0, t) = g^(x)a;-"'(^~^')*. (47) 

Finally, the full (x,^) —t correlation is introduced by "folding in" the Regge ansatz 
for the (x,t) dependence into the DD concept for the {x,C,) dependence. One then 
defines Regge-type DDs : 

F'?(/3, a, t) = F\P, a, 0) /^-"'(i-/^)* (48) 

where F'^{P, a, 0) is given by the form of Eq. (I42l) . To satisfy the polynomiality rule 
of Eq. ( 1T0|) for the GPD H, is has been shown in Ref. [79] that a D-term has to be 
added. The full (x, ^, t) dependence for the GPD H in the VGG model then reads : 



H'^{x, ^,t)= I dad(36{x -{3- ^a)F'>{/3, a, t) 



+ Oi^-\^\)^D(^,t], (49) 



with F'i{P,a,t) defined by Eq. 
(ii) Parameterization of the D-term 



The t-dependence of the D-term is unknown. The D-term, being odd in x, is not 
at all constrained by the FF sum rules of Eq. ([7]). VGG adopts a factorized form 
with a dipole behavior in t with an adjustable mass scale. 
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(iii) Parameterization of the GPD E 

The parameterization of the GPD corresponding with a nucleon hehcity flip 
process, is less constrained as we don't have the DIS constraint for the x-dependence 
in the forward limit. 

One contribution to the GPD E'^ is however determined through the polynomiality 
condition, which requires that the D-term contribution is cancelled in the 
combination H + E. Therefore it contributes with opposite sign to H and E. 
Similarly to Eq. (H9|) for H"^, E'^ is parameterized in the VGG model by adding a 
double distribution part to the D-teim as : 

E'^{x,^,t) = El^{x,^,t) - e{^-\x\)^D(^^,?j , (50) 

where E^j^ is the double distribution part. 

In the forward limit, the DD part reduces to the function e'^{x) = E'[)jy{x,0,0) 
which is a priori unknown, apart from its first moment : 

dx e%x) = k'^ , (51) 

1 

where and k'^ are the flavor combinations of the nucleon anomalous magnetic 
moments given by : 

K« = + = 1.673, k'^ = + 2 k" = -2.033. (52) 

For the x-dependence of the forward GPD e'^{x), a sum of valence and sea-quark 
parameterization is implemented in VGG, according to Ref. [3l] as : 

e"(x) = A^u,ai{x) + B^6{x), 

e\x) = AU,ai{x) + BU{x), 

e'{x) = , (53) 

where the parameters A", A'^ are related to J", J'^ through the total angular 
momentum sum rule, which yields : 



^ - Ml""' ' ^ ^ 

and where the parameters 5", B'^ follow from the first moment sum rule Eq. flSTl) 
as : 

= k" - 2 A", B'^ = k'^ ~ A'^. (55) 

Such parameterization allows to use the total angular momenta carried by u- and 
d-quarks, J" and J"^, directly as GPD fit parameters, and can be used to see the 
sensitivity of hard electroproduction observables on J" and J*^, as will be shown 
further on in section |H 

Starting from the model for the forward distribution e'^(x), the ^-dependence of the 
GPD E1j^{x,C,,0) is generated through a double distribution K'^{f3,a,t) as : 

/■I r^-m 

El^{x,i,t) = df3 da 6{x-f3-aOK%f3,a,t) . (56) 

J-i J-i+Wl 
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The double distribution K'^{f3, a, 0) is taken in analogy as in Eq. (H^ . by multiplying 
the forward distribution e'^(/3) with the same profile function as in Eq. f H2|) as : 

K'i{(],a,t) = h{P,a)e''{P) . (57) 

The parameterization of Eq. ( l53l) yields for the GPD i?^^(x, ^,0) : 

Kd(^,(,0) = -Bl,(i,«,0) 

where the first (second) term originates from the valence (sea) contribution to e'^ 
respectively in Eq. (J53ll . and where the parameter b is the power entering the profile 
function. 

For the t-dependence of the GPD E, the VGG model uses a Regge ansatz, which 
was constrained in Ref. [82] to provide a fit to the Pauli form factor F2. Since 
the large-t behavior of F2{t) goes steeper than the Drell-Yan-West relation 
implies a different large-x behavior of e''{x) compared with q{x). To produce a 
faster decrease with t, a simple ansatz is to multiply qvai{x) by an additional factor 
of the type (1 — x)^«, thus modifying the x ~ 1 limit which is the region driving 
the large-t behavior of F2 (t) , as we discussed previously. This yields for the valence 
part : 

KUP,a,t) = hiP,a)N,q,aim^-Pr'r''''. (59) 

where the normalization constant Ng is determined from the anomalous magnetic 
moment, and where the Regge slope a' and the parameter rjq, which determines 
the large-x behavior of the forward GPD e'^{ to be determined from a fit 

to the nucleon Pauli form factor data as. In contrast, the sea-quark cannot been 
constrained by the F2 FF data. It has for simplicity been assumed to have a same 
t-dependence within VGG as for the valence part. 

(iv) Parameterization of the GPD H 

For the (x,(^)-dependence, the GPD H is also based on the DD ansatz with the 
replacement of the unpolarized PDF q{/3) in Eq. ( |42|) by the polarized PDF Ag(/5), 
so as to obtain the appropriate forward limit of Eq. (E]): 

F{/3,a,0) = h{/3,a)Aq{P) (60) 

For the t-dependence, in principle, a Regge ansatz similar to the one used for the 
unpolarized GPDs (Eq. fj47|) ) could be used. However, at this time, given the 
relatively few experimental constraints from DVCS on this GPD, a t-factorized 
ansatz has been kept: 

H{x, ^,t)= [ dad(35{x - l3 - ^a)F{/3, a, t)GAit)/GAiO). (61) 
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(v) Parameterization of the GPD E 

Following the argument of Ref . [HSl [S^ , it is parameterized by the pion exchange in 
the t-channel, which, due to the small pion mass, should be a major contribution: 



with the asymptotic distribution amplitude (pas is given by (t)as{,z) = 3/4 (1 — 2;^), 
and /iyi(^) is the induced pseudo scalar FF of the nucleon. The contribution of 
Eq. (1621) to the E GPD, corresponding to a meson or qq exchange in the t-channel, 
lives only in the — ^ < a; < ^ region and as such, contributes only to the real part 
of the DVCS amplitude. 

In summary, the VGG parameterization is based on very few inputs: 

• a choice for the PDF which drives the forward limit. By choosing a PDF 
parameterization which take into account the evolution equation, a Q^-dependence 
can be introduced in the GPDs. 

• the parameters and hg-, which drive the (x, ^)-dependence and which are set to 1 
by default. 

• the parameters a' and rjq which drive the t-dependence. In Ref. |82j, the fit to the 
proton and neutron FF data yielded a' = 1.105 GeV~^, 77^=1. 713 and 7^^=0.566. 

• the parameters Ju, Jd which control the normalization of the E GPD and which 
are unknown a priori. 

3.1.4. The GK model The GK parameterization of the GPDs has been developed in 
the process of fitting the high-energy (low x) DVMP data and has been published in a 
series of articles [85| [861 [87] . There are numerous data available for DVMP and since 
the same GPDs as for DVCS enter in the DVMP handbag diagram (Fig. [9]), strong 
constraints on the GPD model parameters can be derived, which are not present in 



• In Ref. [85] the DVMP 2-gluon exchange handbag diagram amplitude (Fig. [9]- 
right) was derived for exclusive p° and cj) electroproduction on the proton, taking 
into account some higher-twist corrections (Sudakov suppression and transverse 
momenta of the quark). The authors proposed a DD-based ansatz for the gluon 
GPDs and compared their calculation to the HERA data. 

• In Ref. [86] the 2-quark exchange handbag diagram amplitude (Fig. [9]-left) was 
added and a DD-based ansatz for the quark GPDs and E'^ was proposed, which 
we will describe in the following. 



^d/p ^ 1 ^(3) 

2 




(62) 



VGG. 
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• In Ref. [S7j, exclusive 7r+ electroproduction on the proton was investigated, which 
allowed to derive a parameterization for the H and E GPDs (as well as for the 
transversity GPD Ht-, which will not discuss here). 

Like VGG, the GK model is based on DDs for the (x, ,^)-dependence. In VGG, 
the h exponents in the profile function of Eq. (H2l) are usually taken as 1 but they are 
essentially unconstrained and left as free parameters due to the lack of constraint from 
the DVCS data. In GK, the h parameters are taken as 1 for valence quarks and 2 for 
sea quaks. This values correspond to the asymptotic behavior of quark and gluon DAs 
respectively. 

For the t-dependence, the GK GPD is expressed (at ^ = 0) as its forward limit 
multiplied by an exponential in t with a slope depending on x : 

GPD\x, e = 0, t) = GPD\x, ^ = 0,t = 0) e*P'(^) (63) 

with a Regge-inspired profile functional form: 

Pi{x) = a[\nl/x + bi. (64) 

The label i stands for valence or sea quark flavours, or gluons. Gluon GPDs are in 
principle taken into account in GK. This is a difference with VGG which takes into 
account only quark (valence and sea) GPDs. However, since the present review focuses 
on the valence region and on the leading-twist leading order domain, the gluonic degrees 
of freedom are not included in the following calculations. 
For quark GPDs, Eq. ( 163|) can be rewritten: 

GPD\x, t) = g(x)x-"'* e''-* (65) 

in order to better compare to the VGG ansatz of Eq. (147|) . 

The GK t-dependence is different from the one of the VGG model in that: 

• There is an x-independent term in the exponential (associated with the parameter 

• The x-dependence of the t-slope has an extra (1 — x) factor in VGG (Eq. ( HTIl ). 

The parameters in Eq. ( l63i) and ( l64l) are determined by the analysis of DVMP 
data in the kinematical region ^ < 0.1, > 3 QeV^ W > 4GeV and -t < 0.6 GeV^ 
The data sensitive mostly to the GPD H are available over a wide range, while the 
existing data which have a higher sensitivity to E, H, E are available only in a restricted 
range. Therefore, a Q^-dependence on the GPD H is taken into account through 
the dependence of the PDF used in the DD ansatz [86l EZj (like in VGG) while it is 
neglected for the E, H and E GPDs. It is also ensured that the valence quark GPDs 
are in agreement with the nucleon form factors at small t and that all GPDs satisfy 
positivity bounds [371 [90]. We now detail the parameterization of each GPD. 

(i) Parameterization of the GPD H 
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The forward limit of the GPD H is the usual unpolarized PDF. To allow an analytic 
evaluation of the resulting GPD, PDFs are expanded on a basis of half-integer 
powers of x: 

3 

H\x, ^ = 0, t = 0) = - xf"'^^ x^'^ (66) 

j=0 

where i represents various quark flavours. The Q^-dependent expansion coefficients 
Cij = CijlQ"^) have been obtained from a fit to the CTEQ6M PDFs [88] and are 



summarized in Ref. [22] • The parameters appearing in the profile functions (I63p 
obey linear Regge trajectories: 

ctHi = aHi{0) + a Hit. (67) 

It is assumed that aHsea.it) = c^Hgit) and aHg{0) = 1 + 5^ as seen in the HERA 
experiments. The expression of the GPD H stemming from the expansion of 
Eq. (ESI) is: 

3 

H,{x, e, t) = e'"^' c^^H,,{x, 1 1) . (68) 

j=0 

where integrals Hij are written down in Ref. ^86j- The slopes bni are modeled by: 

2 

bHsea = hug = 2.58 GeV-2 + 0.25 GeV'^ In ^ , (69) 



Sea quark GPDs are further simplified [87] as: 

Tju rrd TTS 

sea -"sea "'s-"sea ' 

with = 1 + 0.68/(1 + 0.52 In gVQo)- (70) 

The fiavor symmetry breaking factor Ks possesses a Q^- dependence fitted from the 
CTEQ6m PDFs. The parameters in the previous equations are determined by the 
HERA p° and data, 
(ii) Parameterization of the GPD E 

The constraints on E come mostly from the Pauli FF data [90J, through the sum 
rules of Eq. ([7]). A DD ansatz is also used. ii^(x, ^ = 0, t = 0) is parameterized with 
a classical PDF functional form: 



r(2 - a,ai + (3t 



EU^, e = 0, t = 0) = fTT^— ^fT^xlW'^^^"""^^ - • 

1(1- avaljl (1 + P^alJ 

where the ratio of F functions ensures the correct normalization of E at t = 0. 



The fits to the nucleon Pauli form factors performed in Ref. [SS] fix the parameters 
specifying E for valence quarks to /3^^i = 4 and /S^g^j = 5.6. The trajectory a^vai 
and slope parameter bEi are assumed equal to the corresponding H parameters. 
The GK model of the gluon and sea E GPDs has been given in Ref. [91] following 
an idea of M. Diehl and W. Kugler [92]. The DD ansatz is used again and the 
forward limits of the gluonic and strange quark GPDs are parameterized as 

E'{p,^ = t = 0) = Nsp-'-'^{l- pf^^ (72) 
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using /3es=7 and the same Regge trajectory as for H. The sea is supposed to be 
flavour-symmetric. The slopes of the residues Bei are estimated as : 

bEs = 0.9 bHg. (73) 

The normalization Ns of is fixed from saturating a positivity bound for a certain 
range of x [9Tj (which does not allow to fix the sign of Ng) : = ±0.155. 

(iii) Parameterization of the GPD H 

The Bliimlein-Bottcher results [93] are taken to describe the forward limit of if* 
IHEIIHT]. Only Hi^^ is modeled and constrained by the HERMES data [93 [95]. ^sea 
is set to zero. In the same spirit as the modeling of GPDs H and E, the forward 
limit HI^y{x,^ = 0,t = 0) is written following a DD ansatz and in an analytical 
expansion, with the following profile function: 

2 

HU^, ^ = 0,t = 0) = r], (1 - xf ' (74) 

i=o 

where i = u,d. The factors rj^ and rja guarantee the correct normalization of the 
first moment of Hl^^ which is known from F and D values and /3-decay constants 
?7„ = 0.926 ± 0.014 and 77^ = —0.341 ± 0.018. The normalization factors Au and 
are defined by: 



A-' = Bil-as^,4) 



~ ,~ ^-(^m , - (2-ajji)(l -a^i 

QO + Cil ^ r Ci2 



,(75) 



E^e = = 0(1^1 < ^ (^) , (76) 



(6 - o;^i)(5 - a^i). 

where B{a,b) is Euler's beta function. The coefficients c can be foimd in Ref. |89j . 
(iv) Parameterization of the GPD E 

The GPD E is also determined only for valence quarks. Its sea part is set to 0. As 
for VGG, its modeling takes into account the pion pole contribution which reads 
El El]: 

Fpjt) ^ /x + e 
4^ 2^ 

where Fp is the pseudoscalar from factor of the nucleon. The pole contribution to 
the pseudoscalar form factor is written as 

Fp{t) = -rn^fn 5 ■ (77) 

t — 

where denotes the mass of the pion and g-^NN — 13.1 is the pion-nucleon coupling 
constant, /,r is the pion decay constant. 
The pion's distribution amplitude $^ is taken as: 

$^(r) = 6r(l - r) [1 + asCf '(2r - 1)] . (78) 

with 02 = 0.22 at the initial scale Ql = A GeV^. The form factor of the pion-nucleon 
vertex F^atat is described by 



~ ^^2 T \'^) 
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with Atv = 0.44 GeV. Such a hadronic FF is not present in the VGG 
parameterization of E. 

A non-pole contribution, which is not present in VGG, is added and modeled in the 
same way as E and E, i.e. a functional form for the forward limit is assumed, 
then skewed with a profile function in a DD ansatz. Flavour independence of the 
Regge trajectory and the slope of the residue are assumed. The forward limit 
reads [871 196]: 

EUx, e = 0, t = 0) = Ar|a;°^W (1 - xf . (80) 

The following values for the various parameters involved are a£;(0) = 0.48, a'^ = 
0.45, = 0.9 GeV-2, = 14.0 and A^J = 4.0. 

3.2. Dual parameterization 

3.2.1. Evolution equations and conformal symmetry By definition, conformal 
transformations change only the scale of the metric of Minkowski space, and in particular 
leave the light cone invariant. The whole conformal group admits a particular subgroup, 
named collinear conformal group, which maps a given light ray onto itself. This is of 
special relevance for hadron structure functions since in the parton model, hadrons are 
viewed as a bunch of partons moving fast along a direction on the light cone. It helps 
classifying fields according to their collinear conformal symmetry properties. For details 
we refer to the review of Ref . [97] . 

Although QCD is not a scale invariant theory (it exhibits a spectrum of massive 
bound states), conformal symmetry is a symmetry of the classical theory when quarks 
are considered as massless. It is thus relevant for renormalization at leading order 
since the counter terms satisfy the symmetry properties of the tree-level (classical) 
Lagrangian. As operators with different quantum numbers (or symmetry properties) 
do not mix under renormalization, conformal symmetry is a powerful tool to separate 
operators at leading order. 

In particular Gegenbauer polynomials C*™^^ parameterize the local conformal 
operators associated to the twist 2 matrix elements used to define PDFs or GPDs. 
They diagonalize the ERBL evolution equations that describe the evolution of GPDs in 
the inner region — < x < f33l [341 |35l |36l |371 [38] where GPDs probe the presence of 
quark- ant iquark pairs in the nucleon. This is the region of interest when representing 
GPDs as an infinite series of t-channel exchange resonances, as in the case of the dual 
model, or alternatively of the Mellin-Barnes representation. Therefore expanding GPDs 
on a series of orthogonal Gegenbauer polynomials Cn^'^ is an appealing starting point to 
parameterize GPDs. 

3.2.2. Partial wave expansion and CFFs The dual parameterization of the GPDs 
is based on a representation of parton distributions as an infinite series of t-channel 
exchanges [98]. For the unpolarized GPDs, one defines electric (H^) and magnetic 
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(if*^) GPD combinations : 

=H{x,i,t) + -^E{x,i,t), (81) 

H^^'\x, e, t) = H{x, t) + e, t), (82) 

which are suitable for a t-channel partial-wave expansion, which read for the singlet 
combinations as [99] : 

-n..o^.||:.ifw<i-i;)(i-i;)ci(|).(i), 
-n...)^.||:.if'(.)<i-i;)(i-i;)c-(|)i./(i), 



odd even 



odd even 



(84) 

where Cn'^{z) are the Gegenbauer polynomials, Pi{z) are Legendre polynomials, and 
Bni{t) are generalized form factors. Note that for H^^^ intermediate states with l^'-^ = 
0"'""'", 2"'""'", ... contribute, whereas for intermediate states with l^*-^ = 2"*""*", 4'*'"'", ... 

contribute. 

As stated before, the expansion onto a basis of Gegenbauer polynomials allows a 
trivial solution of the QCD evolution equations at leading ordei[§|: The Q^-evolution of 
the generalized form factors Bni{t,Q'^) reads |105] : 

In 

BniM') = B^,{t,Ql) {^^$l^Y\ (85) 



lnQ7A2 



with /3o = 11 — 4^/ and 



^n = \{^+ ( \ - 4(^(ri + 1) + 7^)) , (86) 
3 \ n[n + 1) / 

where denotes the Digamma function and 'Je the Euler-Mascheroni constant. 

At fixed X and ^ the series on the rhs of Eqs. fl83| \8^ are divergent: the sums 



i^l^^ and H^^^ have a support — 1 < a; < +1 while each term of the expansions have 
a support < X < +C,- However, these formal series can be recast onto convergent 
Gegenbauer polynomial expansions. For example the electric singlet GPD reads |100j : 

oo 

iff )(x, e, t) = 2(1 - x') M^, t)Cl'\x). (87) 



odd 



§ In fact QCD evolution equations "commute" with the parameterization Eq. (jBSj): the GPD at some 
input scale Qq has the same form as the GPD at another scale Qj ^ feature that is usually absent from 
the double distribution representation with the factorized Ansatz involving the profile function /i.(/3, a) 
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The coefficients An are defined by: 

Ai^f\- 2n + 3 ^ (p+l)(p + 2) ^ 



{n + + 2) - 2p + 3 ^ ' 

odd even 



(88) 



Here Rnp{C) is a polynomial of degree n: 



r(^^ + i)r(|+p)^ ' ' V 2 '2 ' 2 '2 

with 2-^1 the Gauss hypergeometric function. The convergent expression fl87|) has been 
used explicitely for fitting in a truncated form as explained in Sec. 14.2.31 

The procedure to sum the formal series of Eqs. ( |83| |8^ over orbital momentum I 
through analytical continuation was originally outlined in Ref. |100j . We briefly discuss 
this in the following for the function and for simplicity drop the superscript {E). 
For an analogous discussion of H^^^\ as well as the polarized GPDs H and E, we refer 
the reader to Ref. [99] . 

In order to sum the formal series of Eq. (!83|) . a set of generating functions Q2u{x, t) 
{v = 0, 1, ...) are introduced, whose Mellin-Barnes moments yield the generalized form 
factors Bniit) as [100] : 

Bnn+l-2u{t) = [ dxx''Q2u{x,t). (90) 







The functions Q2u{x,t) are forward-like because at leading order (LO), their scale 
dependence is given by the standard DGLAP evolution equation, so that these functions 
behave as usual parton distributions under QCD evolution. Furthermore, the function 
Qo{x,t = 0) is directly related to the parton densities q{x) measured in DIS |100j : 

Qo{x,t = 0) = [q + q] (x) - ^ f ^ [q + q] (z) . (91) 

The usefulness of the dual parameterization originates when expanding the GPD around 
^ = 0. The functions with higher u are more suppressed for small values of ^. An 
expansion with x fixed to the order C,"^" involves only a finite number of functions Q2ii{x, t) 
with /i < z/. 

Within the dual parameterization for the GPDs, the CFFs entering hard exclusive 
observables can be expressed in terms of forward-like functions. For the combination of 
the CFF of Eqs. f|T5| [T6l). corresponding with the electric GPD of Eq. flHTl) . this is given 

by [inniiini]: 

1 

1 



4mj^ IT J X 



^2x/^-x^-l 



(92) 



Hn^ + -^ER^ = -2 [ -N^^\x,t) 

4771 jy J X 
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+ 



1 



- ^Dit) , (93) 
where A^/- = 3 and where the D-form factor is obtained from the D-term of Eq. fl43p as : 

oo 

Dit) = Y,dnit). (94) 

n=l 

Furthermore, in Eqs. (1 ^ 1^ . the function N^^\x,t) is defined as : 

oo 

N^^\x,t) = J2^"'QS\x,t). (95) 

Analogous relations are obtained of the magnetic GPD combination H^^^ of Eq. f l82p in 
terms of a function A^(*^) which is a sum of forward-like functions Q^^^ ■ The amplitude 
within the dual parameterization of GPDs, given by Eqs. ([221 ESI), automatically satisfies 
a dispersion relation, with the subtraction constant given by the D-FF [101] . We will 
discuss the general dispersion relation approach for the DVCS amplitude in more detail 
in Section [331 

The information contained in the LO DVCS amplitude is in one-to-one 
correspondence with the knowledge of the functions N{x, t) and the D-FF D{t), because 
Eq. f[92|) can be inverted |101j . i.e. the functions can be expressed unambiguously 
in terms of integrals over ^ of Hjm and Ejm- This inversion corresponds to the Abel 
transform tomography |102j . for details see Ref. |101[I103] . This equivalence implies that 
the function N{x, t) contains the maximal information about GPDs that is possible to 
obtain from the observables. Another important feature of the expressions (19 2 [ for 
the amplitude is that one can easily single out the contributions to the amplitude coming 
from the forward parton densities. Indeed, the first term in the sum ( l95l) is given by 
the function Qq which is related to the (t-dependent) parton densities by Eq. (!9T|) . The 
advantage of the dual parameterization is that one can clearly separate the contribution 
of the (t-dependent) parton densities from genuine non-forward effects encoded in the 
functions Q25 Q45 • • •• 

3.2.3. Modeling the forward functions A number of phenomenological studies of DVCS 
observables have been made using the dual parameterization. Most prominently, studies 
involving only the forward function Qq have been made. In such a minimal model, 
the X dependence is parameter free as it is completely fixed by the forward parton 
distributions, merely leaving the t dependence of the GPD to be modeled. One typically 
uses a Regge motivated model to correlate the x and t-dependence of the function 
QQ{x,t)^ analogous as it was discussed above for the double distribution model. 
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Such a minimal model was found to overpredict the data at small and intermediate 
xb: Ref. [ini UnSl Hoe] found that DVCS experiments at HERA (HERMES) were 
overpredicted by roughly a factor of 2 (1.5) respectively. At larger values of xb, for 
DVCS experiments at JLab@6GeV, it was shown that the bulk effect of the DVCS 
beam helicity cross section difference can be understood within such a minimal dual 
model |108] . which we show in more detail below. 

To improve on the description, especially at smaller values of xb-, within the 
dual parameterization requires to go beyond the minimal model by keeping more 
generating functions Q2,Q4,---, and extend to the next-to-leading order accuracy. A 
first step to model the functions Q2 and Q4 has been made using a nonlocal chiral 
quark model |107j or by extracting them from comparison with the double distribution 
parameterization [99] . 

For a comparison with data in the valence, we will be using a model for the LO 
forward-like functions Qo^^ and Qg*^'' as |109j : 



Qf\x,t) 



q+{x,t) + ^e+(x,t) 



X dz 
2.L. ^ 



4m^ 



(96) 

^ (x, t) = ^ ^ [g+ (^, t) + e+ (^, t)] (1 + I) , (97) 

For the forward GPD H^{x,0,t) = q+{x,t), we use a Regge type ansatz : 

g+(x,t) = g+(x)x-"'*, (98) 

with a' = 1.105 GeV~^ fixed from the form factor sum rule |82j. For the forward GPD 
-E+(x,0,t) = e+{x,t), we use an ansatz by expressing the magnetic GPD g+ + as a 
Mellin convolution of g+ with a kernel function, modeled as : 

q+{x,t) + e+{x,t) = ^q^{z,t)c(^^y (l-^y , (99) 

where C is a constant, which is to be determined from the second moments of the GPDs. 

The second moment of the forward parton distribution at t = yields the total 
momentum carried by quarks and anti-quarks : 

M| = [ dxxq+{x,0), (100) 



whereas the second moment of the magnetic GPD combination (g+ + e+) at t = yields 
the total quark angular momentum : 

2J9= f dxx[q+{x,0) + e+{x,0)]. (101) 
Jo 

Eqs. fllOOl IIOII) then allow to express the constant C as : 
2J^ Tia + (3 + 3) 

M| r(a + 2)r(/3 + l)' ^ ^ 
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X 

Figure 21. Parameterization for the forward function (q" +e!|i)(a;, t — Q) for Ju = 0.3, 
and for different values of a,/3 : a = 0, (3 — (dashed-dotted curve); a = 1, 
/3 = (dotted curve); a = 2, (3 = —0.5 (dashed curve). The solid curve shows 
the parameterization for the forward function q^{x,t — 0). 



3. 3. Mellin-Barnes parameterization of GPDs 

3.3.1. Partial wave expansion In this section, we will discuss the Mellin-Barnes 
parameterization of GPDs described in the works of |113[ 1114] . For simplicity, we 
do not write the dependence of GPDs on the momentum transfer t. The method is 
based on making a partial-wave expansion of GPDs. It is analogous in spirit to the dual 
model partial wave expansion explained in Sec. I3.2[ although both representations differ 
on the resummation of this expansion. 

In order to recover the Mellin moments of PDFs when taking the forward limit of 
conformal moments of GPDs, one rescales the Gegenbauer polynomials to define the 
polynomials c„(a;, ^): 

for any integer n. Conformal moments of ^ GPD F {F = H, E, H or E) are then 

defined by: 

Fn{0 = f dxCn{x,OF{x,0- (104) 
The polynomiality of GPDs ensures that F„ is a polynomial in ^ of degree at most n + 1. 
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The rescaled polynomials are orthogonaljjj to the polynomials Pn defined by: 



^ y ^2 J r(|)r(3 + 

The conformal partial wave expansion then reads: 

oo 

F{x,0 = 5^(-l)>„(x,0i^n(0- (107) 

n=0 

This is the common basis of the dual and Mellin-Barnes representations. The left- 
hand-side of Eq. (11071) has support x G [— 1,+1] and the right-hand-side has support 
X G [— ^, +^]. Therefore for |^| < 1 this sum has to be divergent and can be understood 
as a formal definition of conformal moments. It can be resummed by means of the 
Sommerfeld - Watson transform |115j : 

Zi Jq smvrj 

where C is a contour in the complex plane enclosing all non-negative integers (which 
are the poles of j 1/ sinvrj with residues {—ly/n). At this stage, it is only assumed 
that the analytic continuations of the functions p„ and of the moments F„ have no 
singularities inside the contour C. Using the residue theorem one relates Eq. (11081) to 
Eq. (I107p . Since the analytic continuation of a function of a discrete variable is not 
unique, this is a non-trivial step. A justification for it is given below. 

Using such analytic continuations pj and Fj , one can deform the integration contour 
C so that all singularities of conformal moments lie left to a straight line parallel to the 
imaginary axis. If the integrand of Eq. (II 101) decreases fast enough at infinitjj^l, one 
obtains: 

1 I'C+ico 1 

F{x,0 = 77: dj-. :p,{x,Om)^ (109) 

where the real c constant is suitably chosenEl- This large-j behavior is another condition 
that should be fulfilled by the analytic continuations pj and Fj. 

The analytic continuation of Pn can be expressed in terms of hypergeometric 2-^1 
and gamma F functions: 

p.(^, = - |x|)^p, + 9{x - o^Q, , (110) 

II Orthogonality is meant in the fohowing sense: 
r+l 

dxCn{,X,C)P7n{x,£,) = (-l)"(^nm, (105) 



where the factor (—1)" is introduced for later convenience, precisely to write Eq. (jlOSp . 

^ The integrand should decrease fast enough to drop the contour at infinity. Mellin moments should 

also have a sub-exponential growth to guarantee the uniqueness of their analytic continuation thanks 

to Carlson's theorem |116| . 

+ c ~ 0.35 is retained for fitting in Ref. jll3) . 
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where: 



X 



V2 •' I \ ^\ IP I f A \ T\ A \ n r, ' 



However the exphcit calculation of the analytic continuation of conformal moments for 
any value of ^ and any GPD model fulfilling the aforementioned conditions is an intricate 
mathematical question. An explicit general procedure is nevertheless described in the 
case 1^1 < 1 in Ref. [TTl] . 



3.3.2. Compton Form Factors in the Mellin-Bames representation To simplify the 
discussion we restrict ourselves to the case of the singlet GPD if+ and its associated 
CFF defined in Eq. fl2^ . Inserting the Mellin-Barnes representation (11091) for the GPD 
H and permuting the integrals over x and j, the CFF 1-L{^, Q^) reads: 

m,Q') = Y^ P^'rfjT^Utan^^')] [C'^+..]HM,^^% (US) 

where we indicated the hard scale explicitly. Furthermore, the coefficients are the 
conformal moments of the hard scattering kernel C"*", and the dots refer to NLO terms 
proportional to a^- Since we will discuss LO results, we only quote the expression for 
CO (for results at NLO, see Ref. [IH]): 

^ g^'r(| + ,) 

^' r(|)r(3 + jr ^ ' 

Evolution of GPDs can be included as well in this formalism, taking the conformal 
moment (C ® E)j of the convolution of hard scattering C and evolution E operators. 



3.3.3. Modeling of GPD conformal moments In the spirit of the dual model, the 
conformal moments of GPDs can be viewed as the result of t-channel exchanges of 
resonances Rj with angular momentum J, taking into account an effective 7*7i?j 
vertex hj, a propagator with an effective Regge pole a{t) and a smooth profile for 
the distribution amplitude of Rj: 

m,t) = 'j2H^^^ j_\t) {i-j^Y ^'^'''"^'''^^^ ^^^^^ 

where Wigner's 5*0(3) functions |117j are denoted dQ ,^ where i/ = or z/ = ±1 depending 
on hadron helicities. They involve Legendre polynomials when = (electric GPD 
combination of Eq. ( ISTj) ) and derivative of Legendre polynomials when |z/| = 1 (magnetic 
GPD combination of Eq. fl82|)l 

Such modeling of conformal moments has been used in Ref. ^113j to fit to 
unpolarized DVCS data at small at LO, NLO (in the MS and CS schemes) and 
NNLO (in the CS scheme). 
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3.3.4- Modeling of the GPD H{^,^,t) within the quark spectator model With a DD 
representation and a t-dependence inspired from a quark spectator model, the following 
functional form for the GPD H is used for fitting: 



where the parameters n, a{t) and p are a priori known. In the valence case these 
parameters are deduced from Regge u and p trajectories and PDF parameterization. In 
the sea case, the parameterization (11161) is requested to reproduce the small xb fits in 
the Mellin-Barnes representation. 

LO dispersion relations (see next section) are implemented by means of a t- 
dependent subtraction constant parameterized with a normalization constant C and 
a mass scale Mc: 



For the valence part, this leaves thus five free parameters to fit data: M (valence), b 
(valence), r (valence) which respectively control the t-dependence, the large x behavior 
and the skewdness effect of the valence part of H, and C and Mc which respectively 
control the normalization and t-dependence of the D-term. One can use an ansatz 
similar to Eq. 11161 for H which introduces three additional parameters M, b and f. We 
come back to this in Section I4.1.4[ 

3.4. Dispersion relation approach to DVCS: general formalism 

As has been discussed in Section [2], the observables entering DVCS are the CFFs, 
which depend on the GPDs. The CFFs correspond with the real and imaginary 
parts of the DVCS amplitudes, as given by Eqs. (jTSl - 1221) . At a fixed value of the 
momentum transfer t and the external virtuality Q^, the analyticity of the virtual 
Compton amplitude in the energy variable u (or equivalently in the variable ^) allows 
to write down general dispersion relations, which relate these real and imaginary parts. 
The different GPD models and parameterizations detailed above are specific ways to 
implement such dispersion relations. In particular, the double distribution, dual, and 
Mellin-Barnes parameterizations all satisfy the general dispersion relation for the twist-2 
DVCS amplitude which we will detail below. 

When making contact with data, the dispersion relation (DR) approach is a very 
useful tool as it allows to put additional constraints on the data by relating different 
observables. For the non-forward virtual Compton scattering (VCS) on a nucleon, 
described by a spacelike incoming photon and a real outgoing photon, the general 
dispersive framework for the 12 independent invariant amplitudes has been developed 
in Ref. [118j . Such a framework has been successfully applied to VCS data in the region 




(116) 




(117) 
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around threshold and into the first nucleon resonance region, and allowed to extract low- 
energy constants, proportional to nucleon polarizabilities from VCS data, see Ref. [119] 
for a review and details. 

When describing DVCS data, one may also use the dispersive techniques to put 
additional constraints on GPD parameterizations of the CFFs. For the twist-2 DVCS 
amplitude, the dispersive framework for the DVCS amplitude has been discussed in a 
number of different works, see Refs. [I20l [ml [Ml [I22l [Ml |I23]. We will discuss here 
the use of DRs for the DVCS amplitude in some detail for the GPD H. Subsequently, 
we will briefly discuss the amplitudes involving the other GPDs. 

We denote the DVCS amplitude which involves the GPD H by A{^,t) : 



dx H+{x, ^, t) 



+ 



X 



C,+ie X + ^ 



IE 



:ii8) 



which depends on the GPD singlet combination if+ defined in Eq. ( I25l) . The real and 
imaginary parts of the amplitude are related to the CFF introduced in Eq. (12^ . 

as : 

A{i^t)= -m,t). (119) 
To write down DRs, we start by introducing the kinematic (energy) variables : 

u = z/' = (120) 

2mt, 2mx 

which allow to define amplitudes which are either even or odd in u. Denoting the DVCS 
amplitude depending on z/ and t by A(z/, t), the unpolarized DVCS amplitude is even in 
u, i.e.: 

A{u,t) = A{-u,t). (121) 

We can then write down a once- subtracted DR for the amplitude A (assuming one 
subtraction is enough to make it convergent) as : 

'2 '"^ du'^ ImA{u',t) 



ReA{p, t) 



A(0,t) + —V 

IT 



^0 



where a subtraction has been made at v 
the elastic threshold. Using Eq 
i as : 



0, and where vq 
we can rewrite Eq, 



(122) 

Q"^ /2m^ corresponds to 
in a DR in the variable 



ReA{i,t) = m 



2 dx ImA{x,t) 



or equivalently 



ReA{^,t) = A{t) -V dxH+{x,x,t) 



X 



(123) 



(124) 



where the subtraction term (at zero energy) is denoted by A(t). One notices that in 
contrast to the convolution integral entering the real part of the CFF in Eq. flll9p , where 
the GPD enters for unequal values of its first and second argument, the integrand in the 
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DR (spectral function) corresponds to the GPD where its first and second arguments 
are equal. Combining Eqs. f lllSp and f ll24p allows to re-express the subtraction term as: 



A{t)=vf dx [H+{x,^,t) - H+{x,x,t)] 
Jo 



which is independent of C,- When formally taking ,^ = in Eq. fll25p and use time 
reversal invariance : H{x, —x,t) = H{x,x,t) to convert from the singlet GPD to 
the GPD H, one arrives at the sum rule: 



A(t) = -2 / dx-[H{x,0,t) - H{x,x,t)] 
J-i X 



(126) 



As pointed out in Ref. |124j . since both H{x, 0,t)/x and H{x, x,t)/x are even functions of 
X, their singularities cannot be regularized by the principle value prescription, and there 
are no indications that the singularities of both functions cancel each other. However, 
it was shown [124] that the validity of the sum rule of Eq. (I126P can be demonstrated 
by decomposing the GPD into a double distribution (Hdd) part and a D-term (Hd) 
part. The 1/x integrals of the H^d parts (plus distributions) do not contribute to the 
sum rule of Eq. (11261) . The D-term parts are proportional to a 5-function in x as: 

rd.-«(^'') 



X Nf J_i z 

^£(flM=i(,)_L fdz^p:^, (127) 

X Nf z{l - z) 

where D{z,t) is the D-term, see Eq. (1431) . Using Eq. (I127p . one then obtains for the 
sum rule of Eq. (I126p : 



Ni i_i 1-z 

We thus observes that the subtraction term entering the DR for the DVCS amplitude 
A is directly proportional to the D-term form factor. It can be obtained from the 
Gegenbauer expansion of the D-term, Eq. ( H3|) . as: 

A(i) = jrj: dM. 

n odd 

In practise, one can evaluate the dispersion integral in Eq. (I124p by the ordinary integral: 
ReAiU) = m - f dx {H+ix,x,t) - H+{^,U)} + ^ 

-i7+(e,e,t)ln(^^^), (129) 

which is easy to implement in a numerically stable way. Note that for the case where 
H{x, x,0) ~ 1/x for X — )■ 0, the singularity cancels out of the integral. 

For the CFF involving the GPD E, one can write down an analogous sum rule as 
for H. In this case, the subtraction function is given by —A{t), as the D-form factor 
drops out in the sum of H + E. 
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Analogous to the unpolarized DVCS amplitude, which is even in u, one can also 
write down a DR for the polarized DVCS amplitude, which involves the GPD H, and 
which is odd in z/ : 



A{u,t) = -A{-u,t). 
Assuming an unsuhstracted DR for the odd amplitude, allows to write : 



2i> 

ReA{u,t) = —V 

71 



°° , ,ImA(u',t) 



^0 



(130) 



(131) 



Denoting the polarized DVCS amplitude, depending on ^ and t by y4(^, t), the DR then 
reads : 



ReA{i,t) 



1 



V I dx ImA{x,t) 



X — C, X + C, 



(132) 



with 



ImA{x,t) = n H^{x,x,t), (133) 

with the polarized singlet GPD defined as in Eq. 

In Figs. [22] and [231 we show a DR evaluation of the DVCS amplitude corresponding 
with the GPDs H and E respectively, within both DD and dual parameterizations. 
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Figure 22. Comparison of ttie real and imaginary parts of the CFF related with 
the GPD H for the proton at i = 0, excluding the D-F¥ subtraction term. The 
DD parameterization for h^ai — ^sea = 1 is shown by the solid red curves. The dual 
parameterization based on the forward function Qo, Eq. (PT|) . is shown by the dashed 
blue curves. 
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Figure 23. Comparison of ttie real and imaginary parts of the CFF related with 
the GPD E for the proton at t = 0, excluding the ZJ-FF subtraction term. The DD 
parameterization for J„ = 0.3, Jd = 0, and for byai — bsea = 1 is shown by the solid 
red curves. The dual parameterization is shown for — 0.3, Jd ^ 0, and for different 
values of a, (3 in the model for the forward function e+ : q = 0, /3 = (dashed-dotted 
black curves); a — 1, (3 — (dotted green curves); a — 2, (3 ^ —0.5 (dashed blue 
curves). 
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4. Data vs Models and Fits 

4-1. Comparison of models to current data 

In this section, we compare the outputs of the models that we describe in the previous 
section with the currently existing DVCS data that we presented in section 12.31 

Ji^.l.l. VGG vs data We begin by exploring a few configurations and options offered 
by the VGG model. In Fig. |24l we show the JLab Hall A unpolarized (top panel) and 
beam-polarized (bottom panel) cross sections data [5B] as a function of 0. The green 
dotted curve in the top panel shows the result of the BH calculation. It is at the origin 
of the characteristic shape of the cross-section which peaks at low and large 0. However, 
it is clearly not sufficient to fully describe the data. There is certainly room for another 
process, i.e. DVCS as we believe. 

The solid red curve in Fig. |2l] shows the result for the BH+DVCS process when 
only the H GPD, with h^ai = ^sea=l (Eq. fH2|) ) and without the D-term, is included. 
The calculation is now rather close to the data but nevertheless it doesn't perfectly 
describe the distribution. In the —t = 0.17 GeV^ bin, it underestimates the low and 
large data while it gives a good agreement around 0=180°. In the larger —t bins, the 
situation is opposite: it gives a good agreement with the low and large data while it 
underestimates the data around 0=180°. Regarding the beam-polarized cross sections 
(bottom panel of Fig. I25l) . we see that this configuration, with only the H GPD and 
without the D-term, provides a relatively good agreement with the data for the three 
lowest t-bins. This observable is therefore largely dominated by the H GPD, which was 
expected (see Eq. (|27|) ). One notes however a disagreement between the data and the 
calculation for the largest t-bin. This might be a shortcoming of the VGG model in the 
H parameterization but this might also be a sign of higher-twist effects turning in as 
—t increases, as these calculations, we recall, have been done at the leading-twist order. 

The dashed red curve in Fig. [21] shows the result of the same configuration but for 
bvai = &sea=3. The cffcct is to decrease the unpolarized cross section by several percents 
at low —t and the beam-polarized cross section by a couple of percent (in absolute value) 
for all —t values. The effect of these parameters is therefore not dramatic. 

The black dotted curves in Fig. [2^ show the inclusion of the D-term in the previous 
calculation. We recall that the D-term contributes to both H and E GPDs, so that this 
calculation contains also an E contribution. The D-term has a significant influence on 
the cross section: it tends to increase the cross section around 0=180° and thus improve 
the agreement with the data but at the same time it reduces the cross section at low 
and large which is actually not particularly desired. The inclusion of the D-term in 
the calculation has no effect on the beam-polarized cross section (the dashed blue and 
black dotted curves are superposed in Fig. [2l]-bottom panel). Indeed, this observable 
is sensitive only to the imaginary part of the DVCS amplitude (see Eq. ( 1271) ) while the 
D-term, which lives uniquely in the — ^ < a; < ^ region, contributes only to the real 
part. 



GPDs in the valence region from DVCS 



57 




Figure 24. Unpolarized (top row) and beam-polarized (bottom row) cross sections 
for the e~p — ^ e~p7 reaction. Tlie solid circles are the data points from JLab/Hall 
A [58] . The dotted green curve is the result of the BH alone calculation. Four different 
configurations of the VGG model are displayed. The solid red curves are the VGG 
calculation with only the H GPD, without the D-term, and with b^ai = ^sea=l- The 
dashed red curves are the same but with b^ai — bgea—i- The dotted black curves 
correspond to this latter calculation with the addition of the Z?-term. The blue 
dashed curve is the VGG calculation with only the H GPD, without the D-term, and 
with hyai = &sea=l at the slightly shifted kinematics xb=0.365 and Q^=2.35 GeV^ 
(compared to a;B=0.36 and Q^=2.3 GeV^ for the other calculations). 



Finally, we wanted to highlight the importance of determining precisely the 
kinematics of the observables, an issue that we brought up in section 12.21 The values of 
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xb and provided by Ref. [5B] are respectively 0.36 and 2.3 GeV^. Since this leaves 
free the third digit, we explore the kinematics corresponding to xs=0.365 and Q^=2.35 
GeV^, i.e. a set of extreme values yielding the rounding xb=0.36 and Q^=2.3 GeV^. 
The dashed blue curve in Fig. [2^ show the result, for the configuration with only the H 
GPD contribution with h^ai = ^sea=l and without the D-term {i.e. directly comparable 
to the solid red curves). The effect is non-negligible, decreasing the unpolarized cross 
section by up to 15% at the lowest —t values and the beam-polarized cross sections by 
even a larger amount. Keeping this potential effect in mind, we continue our studies in 
the following with the published ^^=0.36 and (5^=2.3 GeV^ kinematics. 

We do not display calculations with the E GPD because its only effect on the 
observables of Fig. |24]is to increase by 1 or 2% the unpolarized cross section (it cannot 
contribute to the beam-polarized cross section). 

Having discussed these few effects on the H GPD, we are now going to focus on 
the effect of the H and E GPDs. In the following figures, we will keep h^ai = &sea=l for 
H. We will compare the VGG calculations to the JLab Hall A, CLAS and HERMES 
data. We will show systematically four sets of curves and configurations: 

• Only the GPDs H contribution, without the D-term. In the following figures, this 
configuration will be described by the red solid curves. 

• Adding, with respect to the previous configuration, the H contribution. This 
configuration will be described by the dashed red curves. 

• Adding, with respect to the previous configuration, the E GPD with its valence 
and sea contributions (see Eq. ( l53l) ). The values of (J„, Jd) are taken as (0.3, 0.). 
The D-term contribution to H and E is included. In the figures, this configuration 
will be described by the black dash-dotted curves. 

• Changing, with respect to the previous configuration, the values of (J^, Jd) which 
are taken as (0., 0.3). In the figures, this configuration will be described by the 
blue dashed curves. 

Fig. [25] shows the results of these four calculations for the JLab Hall A unpolarized 
and beam-polarized cross sections. By comparing the dashed red curves [H+H 
contribution) to the solid red curves [H-onlj calculation), one sees that the H GPD 
has very little effect on the unpolarized cross section as the two curves are barely 
distinguishable in the top panel of Fig. [2H The effect of H on the beam-polarized cross 
section is more visible. It tends to increase by ~ 15% the amplitude of the sine-like 
modulation. 

One sees that the VGG E GPD has almost no influence on the beam polarized 
cross section as the black dash-dotted curves and the blue dashed curves are essentially 
superposed on the dashed red curves in Fig. [251-bottom panel. The beam polarized cross 
section is thus largely dominated by H with a small extra contribution of H. Concerning 
the unpolarized cross section, the sizable influence of the GPD E is essentially through 
the D-term as the dashed-dotted black curves and dashed blue curves of Fig. [25]- 
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top panel are almost indistinguishable. In other words, there is no sensitivity of the 
unpolarized cross section to the ( J„, Jd) contribution of the E GPD. 




Figure 25. Unpolarized (top row) and beam-polarized (bottom row) cross sections 
for the e^p e~pj reaction. The solid circles are the data points from JLab/Hall 
A [58]. Four different configurations of the VGG model are displayed. The solid red 
curves are the VGG calculation with only H. The dashed red curves are the same 
calculation with the addition of the H contribution. The dashed-dotted black curves 
are the same with the addition of the E GPD with its valence and sea contributions 
with (J„, Jd)=(0.3,0.) and the D-term. The dashed blue curves are the same but with 
( Ju, Jd)=(0.,0.3). The dotted green curve is the result of the BH alone calculation. 

In Fig. [261 W6 compare the four VGG calculations to the beam spin asymmetries 
of the CLAS collaboration. We recall that these asymmetries are the ratio of the beam- 
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polarized to the unpolarized cross sections. Since the VGG calculation with only H 
and with or without H is in general underestimating the unpolarized cross section while 
describing correctly the beam-polarized cross section at low — t, as we saw in Fig. [251 it 
should be expected that the beam spin asymmetry be overestimated at low —t. This is 
indeed what we observe in Fig. |26l Since the addition of H increases the amplitude of the 
beam polarized cross section (Fig. [251-bottom panel), the corresponding asymmetry is 
also amplified. At larger —t (>^ 0.8 GeV^), where the leading-twist handbag formalism 
is expected to be less valid, the agreement between the data and the calculation is 
better. However, we saw in Fig. [25] that, for the largest —t bin, the VGG calculation 
was overestimating the beam-polarized cross section so that the better agreement for the 
beam spin asymmetry might simply result of the ratio of two overestimated quantities. 
This clearly shows the limit of comparing calculations to one single asymmetry. 

Adding E in the VGG calculation (black dash-dotted and blue dashed curves in 
Fig. [26l) moves the VGG beam-spin asymmetries calculations closer to the data. This 
could be expected since the addition of E tends to increase the cross section around 
(/)=90° (see Fig. [251) . However, some discrepancy clearly remains. One can also note 
that in general, as increases, the agreement between the calculations and the data 
tends to improve. 

We finally compare in Fig. [271 our four VGG calculations to the lower xb HERMES 
domain. We show in this figure the nine asymmetry moments which are expected to 
be non-null in the leading-twist handbag formalism. We added as a tenth observable 
the A^^j^'^ moment (bottom right plot in Fig. [27l) which is expected to be small 
(suppressed by powers of Q) in this approximation. However the data shows a rather 
large asymmetry which cannot be described by any leading-twist handbag calculation, 
whatever the parameterization of the GPDs. 

In Fig. [271 we see that the main trends of the data are reproduced by the VGG 
model: for instance, as —t increases, the trend towards increasing negative values of 
Ac and the rise of A^^^'^ and Auyj, etc... Except for the amplitudes of Aq and 

A^^^"^ which are overestimated, the VGG calculation provides a good description of the 
amplitudes of the nine leading-twist observables. We note the particular sensitivity of 
Ac and A™'^'^ to the E GPD. It actually mostly comes from the D-teim contribution to 
E, since there is httle difference between the dashed-dotted black and dashed blue curves 
(non-D-term contribution to E). However, since Ac and A'^^^'^ are largely overestimated, 
it seems that no really reliable conclusion on E or the D-term can be extracted at the 
moment. 

The transversally polarized target asymmetries are also expected to be particularly 
sensitive to E and one does see some difference between the two ( J^, J^) configurations. 
However, we see that these observables are actually largely dominated by H and that 
E comes only as a small variation around H. It is also difficult in those conditions to 
extract a reliable information on E, as long as H is not determined a the few percent 
accuracy. 

As expected, A^']^'^ is particularly sensitive to the H GPD contribution, which is 
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Figure 26. Beam spin asymmetries at (f)=90° as a function of t as measured by the 
CLAS collaboration [61] with VGG calculations. The convention for the curves is the 
same as in Fig. [25] 



necessary in order to explain the magnitude of the data. Finally, the large amplitude of 
the A^^j^*^ moment is a puzzle. 

4.1.2. The GK model vs data We have ran the GK model in three different 
configurations: 

• Keeping only the GPD H. In the figures, this configuration will be described by 
the red sohd curves. 

• Adding with respect to the previous configuration the contribution of the H GPD. 
This configuration will be described by the dashed red curves 
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HERMES data with VGG mode 
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Figure 27. Ten azimuthal moments as a function of —t as measured by the HERMES 
collaboration [MllMllMllSaiMllMllZQlin] with VGG calculations. The convention 
for the curves is the same as in Fig. [25l 



• The "full model", i.e. with the contributions of all four GPDs. This configuration 
will be described by the black dash-dotted curves. 

Fig. [281 shows the results of these three calculations for the unpolarized and beam- 
polarized cross sections as measured by the JLab Hall A data [58] . 

We observe some features very similar to the VGG calculation. The beam-polarized 
cross section (bottom panel of Fig. [28]) is well described for the three lowest —t bins 
and is mostly the result of the H GPD contribution. The inclusion in the calculation of 
the H GPD increases, like for VGG and in about the same proportions, the amplitude 
of the polarized cross section. Since the dashed-dotted black curves ( "full model" ) are 
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essentially superposed on the dashed red curves, we conclude that E (and E) have 
basically no influence on this observable. In the largest —t bin, the GK calculation 
tends to overestimate the data although to a lesser extent than VGG. 

For the unpolarized cross section (top panel of Fig. |28|) . we also see features similar 
to the VGG model. The calculation with only the H GPD is rather close to the data but 
is not fully satisfactory. There is clearly some missing strength for (j) around 180°, where 
the handbag DVCS contribution is expected to be dominant. Adding the contribution 
of the other GPDs barely makes a difference. It might even slightly decrease the cross 
section around 0=180°, worsening the situation. We recall that the GK model has no 
D-term implemented. 

Fig. [29] compares the CLAS beam spin asymmetries with the GK calculations. As 
could be anticipated, the beam spin asymmetries are in general overestimated at low 
—t since the unpolarized cross sections are underestimated (top panel of Fig. I29l) and 
the beam-polarized cross sections are correctly reproduced (bottom panel of Fig. |29|) . 
Adding the contribution of the GPDs other than H tends to increase the disagreement 
between the calculation and the data. This can be attributed to the decrease of the 
unpolarized cross section that we noted in Fig. [281-top panel. However, like for VGG, 
we observe that the GK calculation provides a better agreement with the data for the 
largest values, in which case, the inclusion of the GPDs other than H does help. 

Fig. ED] compares the HERMES azimuthal moments with the GK calculation. Like 
for VGG, the general trend of the data is correctly reproduced. A notable difference 
is that the amplitudes of the Ac and A^^^'^ moments agree much better with the data. 
This can in part be attributed to the absence of a D-term in the GK calculation as we 
saw that; in the VGG model, it was a major contribution, both to H and E. We notice 
again that iheAuyj moments are largely dominated by the contribution of the H GPD 
and that any reliable extraction of the comparatively small E contribution requires a 
control of H at the few percent level. Like in VGG, the strong A^^^^"^ moment is a 
mystery. 

4.. 1.3. The dual model vs data We have ran the dual model in three different 
configurations: 

• Keeping only the GPDs H, without D-term, i.e. for di = 0. In the figures, this 
configuration will be described by the red solid curves. 

• Adding, with respect to the previous configuration, the E GPD with di = —4/3 
and {Ju, Jd)= (0.3, 0.). The {a,f3) values in the model for the function e+(a;,t) 
(Eq. ( l99l) ) are taken as (0., 0.). In the figures, this configuration will be described 
by the black dash-dotted curves. 

• Changing, with respect to the previous configuration, the values of (a, (3) (Eq. ( l99l) ) 
to (2., -0.5). In the figures, this configuration will be described by the blue dashed 
curves. 
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JLab Hall A data with GK model 
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Figure 28. Unpolarized (top row) and beam-polarized (bottom row) cross sections 
for the e~p — > e^p7 reaction. Tlie solid circles are the data points from JLab/Hall 
A [58]. Three different configurations of the GK model are displayed. The solid red 
curves are the GK calculation with only the H GPD. The dashed red curves are the 
GK calculation with in addition the H GPD. The dashed-dotted black curves are the 
"full" GK calculation, i.e. with the contribution of all GPDs. The dotted green curve 
is the result of the BH alone calculation. 

Fig. [31] shows the results of these three calculations for the unpolarized and beam- 
polarized cross sections as measured by the JLab Hall A data [58] . 

The beam-polarized cross sections (bottom panel of Fig. [31]) . are very well described 
for the three lowest —t bins. This gives support to the description of the CFF in 
the dual model. One can note that there is a little influence (a few percents) of the E 
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JLab CLAS data with GK mode 
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Figure 29. Beam spin asymmetries at (/)=90° as a function of t as measured by the 
CLAS collaboration [61 with GK calculations. The convention for the curves is the 
same as in Fig. [28l 



GPD on this observable. Like for the previous models (VGG and GK), the largest —t 
bin shows a discrepancy between the calculation and the data. 

Regarding the unpolarized cross section (Fig. [SU-top), we see that having only the 
GPD H for the DVCS contribution is not sufficient to describe the data. We note that 
the GPD H in the dual model brings less strength to the unpolarized cross section 
than in the VGG and GK models (see Figs. [25] and [281) . Like for the VGG and GK 
models, adding a E GPD contribution doesn't improve the situation: it increases the 
cross section around </)=180° but decreases it at low and large 
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HERMES data with GK r-ode 
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Figure 30. Ten azimuthal moments as a function of —t as measured by the HERMES 
collaboration [BllESlIMllSIllMlIMllIQlEl] with GK calculations. The convention for 
the curves is the same as in Fig. [28l 



In Fig. [321 we compare the three dual calculations to the beam spin asymmetries 
of the CLAS collaboration. We observe the same features as with the VGG and GK 
calculations, i.e. the best agreement between the data and the calculation for the couple 
of bins with > 3 GeV^ bins but otherwise an overestimation (by 20 to 30%) of the 
data. This is obviously due to the underestimation of the unpolarized cross section. 
Adding a contribution of the E GPD can bring the calculation a bit closer (or further) 
to the data according to the parameterization chosen but this is clearly not a main 
factor. 
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JLab Hall A data with dual model 
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Figure 31. Unpolarized (top row) and beam-polarized (bottom row) cross sections 
for the e~p — > e^p^ reaction. The sohd circles are the data points from JLab/Hall 
A [SS]. Three different configurations of the dual model are displayed. The solid red 
curves are the dual parameterization for the vector amplitude (GPD H) only (without 
D-term, i.e. for di =0). The other curves are the dual parameterization including 
both GPD's H and E, for di = —4/3, J„ = 0.3, Jd ~ 0, and for different values of a, (3 
in the model for the function e+{x^t) : a = 0, /? = (dashed-dotted black curves), 
and a = 2, /3 = —0.5 (dashed blue curves). The dotted green curve is the result of the 
BH alone calculation. 
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JLab CLAS data with dual mode 
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Figure 32. Beam spin asymmetries at (/)=90° as a function of t as measured by the 
CLAS collaboration [HI] with the dual model calculations. The convention for the 
curves is the same as in Fig. [211 



In Fig. [331 we compare the three dual calculations to the azimuthal moments 
measured by the HERMES collaboration. We observe the same trends as for the VGG 
and GK models, i.e. a general good description of the trend of the data. However, 
in the details, the Ac and A'^^'^'^ amplitudes are overestimated. The strong sensitivity 
of the A'^f^y^ moment of the E GPD is again noticed. Here, the data clearly favor the 
a = 2, P = —0.5 configuration for the parameterization of the e^{x,t) function. 
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HERMES data with dual mode 
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Figure 33. Ten azimuthal moments as a function of —t as measured by the HERMES 
collaboration [BllESllMllSIlEHlIMllZniIIl] with the dual model calculations. The 
convention for the curves is the same as in Fig. 1311 



4-1 .4- The KM model vs data The KM model is the model originally developped by K. 
Kumericki and D. Mueller |113j based on the Mellin-Barnes parameterization of GPDs 
which we discussed in section 13.31 The parameters of the KM model are determined 
from fitting Hl/ZEUS, HERMES (with or without target polarization data) and JLab 
data (CLAS and with or without Hall A). In the following, we have used three versions 
of the code: 

• KMIO: in this version, the four GPDs are considered: H and H are modeled for 
their valence part by Eq.( I116p . the E contribution is only through the D-term 
and E by the pion pole (like in VGG and GK). There are fifteen free parameters. 
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They are h (val), r (val) and M (val) for the valence part of H (see Eq.( I116p ). 
three similar ones h (val), f (val) and M (val) for the valence part of H, C and Mc 
for the D-term (see Eq.( 11171) ). two for the pion pole (one for the normalization 
and one for the t-dependence) and five for the small x behavior of H: M (sea) 
which controls the t-dependence trough a dipole ansatz, S2 (sea), S2 (gluon) , S4 
(sea) and S4 (gluon) which control the normalization and evolution flow of the sea 
and gluon contributions. One should however note that the five latter parameters 
are determined by the collider experiments which are then just propagated to the 
valence region. In other words, the JLab and HERMES data that we discuss in this 
review only impact the first ten parameters and the KMIO model can effectively be 
considered as having 10 free parameters in this framework. These parameters are 
determined by fitting the JLab Hall A, CLAS beam spin asymmetries and HERMES 
data, excluding the polarized target data. The values of the free parameters and 
more details in the ingredients of the model are given explicitely in Ref. |lllj . In 
the following figures, this version will be described by the black dot-dashed curves. 

• KMlOa: compared to the KMIO version, this version sets H to zero and fixes 
the pion pole, which removes five free parameters and therefore reduces the total 
number of free parameters to five. These are determined by fitting only the 
HERMES (without the polarized target data) and CLAS data, i.e. excluding the 
JLab Hall A data of the fit. The JLab Hall A cross sections are indeed notoriously 
difficult to describe, as was illustrated by the comparison of the VGG, GK and 
dual models with data shown in the previous sections. We will see shortly that 
this "minimal" version of the KM model is indeed unable to describe or predict the 
JLab Hall A cross sections. In the figures, this version will be described by the red 
solid curves. 

• KMM12: this latest version of the KM model includes the HERMES polarized 
target data in the fit (as well as the JLab Hall A data) and has the same ten free 
parameters than KMIO. The values of these parameters and the detailed ingredients 
of this version of the model are given in Ref. |112] . In the figures, this version will 
be described by the blue dashed curves. 

We mention that the code for the KMIO and KMlOa versions of the model can be 
found on the public link of Ref. |110] . We are thankful to K. Kumericki and D. Mueller 
for having provided us with a private version of the KMM12 code. 

The results of the three KM model versions are displayed in Fig. [3l]in comparison 
to the JLab Hall A unpolarized (top panel) and beam-polarized (bottom panel) cross 
sections data. The solid red curves (KMlOa model), which don't include the JLab 
Hall A data in the determination of the parameters are essentially the result of the H 
contribution (with some E contribution through the D-term and E contribution through 
the pion pole, as discussed earlier). We then find results similar to those obtained with 
the three previous models that we discussed, i.e. that the unpolarized cross sections 
are underestimated while the beam-polarized cross section are approximatively well 
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described. It is only with the addition of a if contribution (KMIO or KMM12 models) 
that an agreement is obtained for the description of the unpolarized cross section. An 
issue is that this H contribution is about a factor 3 larger than values given by standard 
parameterizations, such as in VGG or GK. In the KM model, H could therefore be 
viewed as an effective GPD contribution, not clearly linked to polarized PDFs. The 
inclusion of the HERMES polarized target data in the fit (blue dashed curve) does not 
strongly affect the description of the Hall A data (one essentially notes a change of ~ 
10% for the unpolarized cross section lowest | t \ bin). 

In Fig. |35l we compare the results of the KM model versions to the CLAS beam spin 
asymmetries. The comparison between the data and the calculations is available only 
on some limited range, i.e. > 1.5 GeV^ and —t << Q^, due to the restrictions in the 
KM code. We note that those are actually very reasonable limits for the application of 
the GPD formalism in DVCS and should probably be valid for all models, not only KM. 
All three versions of the KM model describe relatively well the data which is expected 
since these data are included in the fit of the parameters for all three configurations. 
The behavior of the three model versions differ somewhat but the data do not allow to 
favor one more than the others. 

In Fig. |36l we finally compare the results of the KM model versions to the HERMES 
azimuthal moments. The polarized target data are described only by the KMM12 
version of the code, which explicitely took those data in the fit. It achieves a relatively 
good description of the nine leading-twist asymmetry moments. As usual, the A^^^^'^ 
is not explained but since it is not related to the DVCS leading-twist formalism, this 
should not come as a surprise. 

In summary, the KM model, in particular the KMM12 version, is, to this day, the 
only model available on the market which achieves a relatively good description (with an 
overall normalized of ~ 1-55 for 95 data points |112j ) of all presently available DVCS 
data in the valence region (and even beyond, at lower xb values). The challenge lied 
particularly in the description of the JLab Hall A data, for which the other three models 
that we discussed failed to give a satisfying description. Though, the price to pay has 
been to introduce in KM a strong H contribution, which remains to be understood. 

4.2. OFF fits 

In the previous section, we have compared the JLab Hall A, CLAS and HERMES DVCS 
data with four models with adjustable parameters: VGG, GK, the dual model and KM. 
Although many general trends of the data are reproduced by these four models, none 
can claim to have a perfect global description of all data with fully reliable inputs (in 
particular, the meaning of the strong H contribution in KM, which achieves the best 
fit all the data sets, remains to be understood). Specifically, the observable which 
is the most challenging to reproduce is the JLab Hall A unpolarized cross section. The 
question arises if this deficiency is due to the limitations of the models which imposes 
some particular functional form for the (x, ^, t)-dependence of the GPDs that might be 
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Figure 34. Unpolarized (top row) and beam-polarized (bottom row) cross sections 
for the e^p — > e^p7 reaction. The dotted green curve is the result of the BH alone 
calculation. The solid circles are the data points from JLab/Hall A [5S]. Three different 
versions of the KM code are displayed. The red solid curves are the results of the 5-free 
parameter model version with the Hall A data and the HERMES polarized target data 
excluded from the fit (version KMlOa). The black dot-dashed curves are the results 
of the 10-free parameter model version with the Hall A data included in the fit, but 
not the HERMES polarized target data (version KMIO). The blue dashed curves are 
the results of the 10-free parameter model version with the HERMES polarized target 
data included in the fit (version KMM12). 
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JLab CLAS data with KM mode 
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Figure 35. Beam spin asymmetries at (/)=90° as a function of i as measured by the 
CLAS collaboration [HT] with the KM model calculations. The convention for the 
curves is the same as in Fig. [Ml 



too constraining or if the data simply don't lend themselves in general to a leading-twist 
and leading-order handbag formalism. 

To give some element of response to this question, we now present an alternative 
way to work on the data. Instead of starting from a model whose parameters, in the 
frame of a particular functional form, are to be fitted to the data, another approach is to 
directly take the CFFs as free parameters and fit them to the data for each individual 
kinematic point. As we saw with (see Eqs.( [271130]) for a few examples), there are 
well established relations (at least at the leading-twist and leading order) between the 
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HERMES data with KM mod( 
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Figure 36. Ten azimuthal moments as a function of —t as measured by the HERMES 
collaboration with the KM model calculations. The convention for the curves is the 
same as in Fig. |34l 



CFFs and the observables and they depend only on kinematical factors. At a given 
{Ee,XB,Q'^,t,(f)) experimental point, these kinematical factors can be determined and 
the idea is thus, for each {Ee,XB,Q^,t,4)) experimental point, to fit simultaneously the 
various observables available for this particular kinematics, taking the CFFs as free 
parameters. The main advantage of this approach is to be model-independent (in the 
limit that the leading-twist and leading order relation between the observables and the 
CFFs are valid) , as the CFFs can vary freely without being constrained or restrained by 
a particular function. The main shortcoming of this approach is that CFFs are fitted, 
i.e. not GPDs themselves. We recall that CFFs, which depend only on C, and t while 
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GPDs depend on x, ^ and t, are (weighted) integrals of the GPDs (over x) or GPDs at 
the hnes x = Therefore, in order to access GPDs, there is, in a second stage, an 
unavoidable additional model-dependent deconvolution to carry out. Nevertheless, we 
will see in the following that already extracting CFFs, in a model-independent way as 
much as possible, as a first step towards GPDs, carries a lot of interest. 

We describe in the following the work in this direction of Refs. |126l 11271 11281 1129[ 
[T251ITT2] . 

4-2.1. "Brute force" least square minimization The method was pioneered in 2008 in 
Ref. |126] . Knowing the well-established BH and DVCS amplitudes (see Eq.(|T3l) for the 
latter process), which provides the relation between the observables and the CFFs, the 
procedure consists in fitting, at each {xb, Q^, —t) kinematic point, the distribution 
(or the moment) of all the available observables at this kinematic point, taking the eight 
CFFs as free parameters. In Refs. \12Q\ \127\ \12S\ 1129] , actually only seven CFFs were 
considered: Hrs, E^e, Hjie, Ejie, Hjm, Ejm and Him- In this work, the eighth CFF 
Ejrn has been set to 0. The reason is that, as was seen in section [3l it is common to 
model the E GPD by the pion pole, which contributes only to the real part of the DVCS 
amplitude. This is essentially the only model assumption in this procedure. Otherwise, 
the other CFFs are free to vary within a 7-dimensional hypervolume, which is only 
bounded by conservative limits: ±5 times the values of the VGG CFFs. We recall that 
some GPDs have to satisfy a certain number of normalization constraints. These are all 
fulfilled by the VGG model and it should be clear that ±5 times the VGG CFFs make 
up conservative bounds. 

It is clear that fitting only one observable with seven free parameters does not 
converge. All data can be fitted with high quality but many combinations of the seven 
CFFs provide an equally good fit and no information can really be extracted on any 
CFFs. However, it was striking to observe in Ref. |126] that fitting simultaneously two 
observables, namely the unpolarized and the beam-polarized cross sections of Hall A, 
with the seven CFFs as free parameters, resulted in a convergence for two CFFs, i-e. 
Hfie and Hj^- This resulted for the first time in quasi-model-independent constraints 
on the Hfie and Him CFFs for the Hall A kinematics. Fig. |37] shows the resulting fits 
of the Hall A data and Fig. 138]- left column the resulting Huf. and Him CFFs obtained, 
displayed clS db function of —t. 

The reason for the convergence of the particular Hue and Him CFFs is, as we saw 
in section 14.11 that the unpolarized and the beam-polarized cross sections are largely 
dominated by these two CFFs (respectively). The fitting procedure allows to determine 
central values for these two CFFs which are the values which minimize the between 
the theory and the data, and two error bars, which correspond to + 1- These error 
bars are asymmetric which reflects the non-linearity and undersconstrained nature of 
the problem (we recall that CFFs enter as bilinear combinations in a cross section). 
The error bars that are obtained reflect actually not the statistical accuracy of the 
data (which are precise at the few percent level) but the influence of the five other 
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CFFs which are sub dominant and do not converge to a particularly well defined value 
within the 7-dimensional hypervolume defined previously. The error bars are therefore 
correlation error bars. If, guided by some theoretical considerations, one can reduce the 
number of CFFs entering the fit as free parameters (like it is done for Ei^n or like it 
will be done in the next section by keeping only the H GPD) or if one can reduce their 
range of variation into an hyperspace smaller than ± 5 times the VGG CFFs as has 
been done here, the error bars can obviously only diminish. In this case, one has clearly 
to make sure that the assumptions are well founded, otherwise, errors will of course 
be underestimated. The present approach should be considered as a most conservative 
estimation of uncertainties and as minimally theory-biased. 

In these conditions, with this fitting algorithm, it was possible to determine : 

• as we just discussed, the Hj^ and Hr^ CFFs at < xb >~ 0.36, and for three 
t values, by fitting simultaneously |126] the JLab Hall A proton DVCS beam- 
polarized and unpolarized cross sections |58] (see Fig. [38]- left column). For the 
lowest —t values, the fitting procedure could not identify a central value for H^e with 
well-defined error bars and so we display only Him- One should also notice that the 
fit of the largest — t-bin is not perfect (x^ ~ 3) as can be seen in Fig. [371 Therefore, 
the small error bar obtained for Hjie at — t =0.33 GeV^ might be underestimated, 
the meaning of an error bar on a fitting parameter for a bad fit being not 
straightforwardly interpretable. 

• the Him and Him CFFs, at < xb >~ 0.35 and < xb >~ 0.25, and for three t 
values, by fitting simultaneously |128j the JLab CLAS proton DVCS beam-polarized 
and longitudinally polarized target spin asymmetries [61], [63] (see Fig. [38]-center 
column). We recall that Acr^iy, and consequently Am is dominated by the Him 
CFF (Eq. ([27])) and that Aajy^,, and consequently Aul, is dominated by the Him 
CFF (Eq. ([281)). 

• the Him, Hue and Him CFFs, at < >~ 0.09, and for four t values, by fitting 
simultaneously |127l 1129] the series of HERMES beam-charge, beam-polarized, 
transversely and longitudinally polarized target spin asymmetry moments [671 [70] 
[66] [68] (see Fig. [38]-right column). In a nutshell, Ac constrains Hrs, Am constrains 
Him and Au^ constrains Him- Unfortunately, in spite of the quasi-complete set of 
observables measured by HERMES, due to unsufficient precision in the data, this 
approach didn't allow to constrain the other CFFs (while in principle, with "ideal" 
infinitesimal resolution, they should). We recall that we showed in Section [^TTl that, 
for instance, the E GPD was actually entering only as a relatively small variation 
with respect to the H contribution in most observables and that without a precise 
determination of H, it is not surprising that no significant information on E can 
be extracted. 

In Fig. [38] the results of all these fits are compiled and shown as empty squares 
(along with and model curves and the rsults of the other fitting strategy that we discuss 
in the remaining of this subsection) . Although error bars are large and the kinematics are 
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E^=5.75 GeV, Xg=0.36, Q^=2.3 GeV^ 
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Figure 37. Result of the fit of the unpolarized (top panel) and beam-polarized 
(bottom panel) cross sections of the e^p e^pj reaction by the fitter code of Ref. |126] 
leaving seven CFFs free. 

not exactly the same between CLAS and the Hall A, it is interesting to note in the case of 
Hjm, that one extracts compatible values from fitting different observables (unpolarized 
and beam-polarized cross sections for Hall A and beam spin and longitudinally polarized 
target asymmetries for CLAS). One can also note that in this method a better precision 
on the extraction of Hjm is achieved by fitting two asymmetries than two cross sections. 
In Fig. [381 some general features and trends can be distinguished: 

• Concerning Hjm, it appears that, at fixed —t, it increases as xb decreases {i.e. 
going from JLab to HERMES kinematics). It is actually possible to extract Him 
at the quasi-common value of —t ^ 0.28 GeV^ from the JLab Hall A, CLAS and 
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Figure 38. The iJ/m, -ffi?,e and Him CFFs as a function of —t for tliree different 
xb values. The empty squares show the results of the 7 CFFs free parameters fit of 
Refs. [126l[T27l[T28l[Tl9]. The solid circles in the HERMES column show the resuh of 
the linear mapping fit discussed in section 14.2.21 (Ref. '112' ) . The solid squares in the 
JLab Hall A column show the result of the 7J-only CFF fit discussed in section 14.2.31 
(Ref. |125| ). The solid red curves show the result of the VGG model (with hyai = ^sea=l 
and without any D-term for H). The solid dashed blue curves show the results of the 
model-based fit of Ref. [113) . 



HERMES data (with a slight interpolation in some cases). We then see in Fig. ^\ 
the xs-dependence of if/^ at —t ~ 0.28 GeV^. One observes the rise of this CFF 
as decreases which is similar to the rise of PDFs (due to sea quarks). We recall 
that Him reduces to a PDF at ^=0 and t = GeV^. In this figure, we also display 
the prediction from the VGG (solid red curve) and KM (dashed blue curve) models 
for Him CFF at -t = 0.28 GeV^. 

• Another feature concerning Him is that its t— slope seems to increase with xb 
decreasing. We recall that the t-slope of the GPD is related to the transverse 
spatial densities of quarks in the nucleon via a Fourier transform (see Eq. (jH])). This 
evolution of the t-slope with xb could then suggest that low- a; quarks (the "sea") 
would extend to the periphery of the nucleon while the high-x (the "valence" ) would 
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tend to remain in the center of the nucleon. We will come back to this discussion 
in section [6l 

• Hjie has a very different t-dependence than if/^ both at JLab and at HERMES 
energies: while Him decreases with —t increasing, Hr^ increases (at least up to —t ~ 
0.3 GeV^) and may even change sign, starting negative at small —t and reaching 
positive values at larger —t. The VGG model (red solid curve), as well as the 
results of the other fitting strategy that we will discuss in the next subsection (solid 
squares) show or tend to show this "zero-crossing" at JLab kinematics. However, 
the KM model does definitely not. Given the large error bars on the fitted Hue-, one 
cannot at this stage clearly favor or exclude any of the VGG or KM models. These 
two models have drastically different predictions for this CFF and more precise 
data on Hrc gerly asked for. 

• Concerning Him-, we notice that it is in general smaller than Him-, which can be 
expected for a polarized quantity compared to an unpolarized one. There is very 
little xb dependence The t-dependence of Him is also rather flat. The weaker t- 
dependence of Him compared to Him suggests that the axial charge (to which the H 
GPD is related) has a narrower distribution in the nucleon than the electromagnetic 
charge. We remark that the slope of the axial FF (the first x-moment of the Him 
GPD) is also well known to be flatter compared to those of the electromagnetic 
FFs. It is very comforting that by studying two relatively different experimental 
processes (DVCS and for instance vr production from which the axial FF is in 
general extracted), one finds similar features. One can also note that there is very 
little Xb dependence for Him ■ 

4-2.2. Mapping and linearization In Ref. [112J, a more elegant method has been 
developped. It consists in establishing a set of relations associating the DVCS 
observables to the CFFs. This is called "mapping". Given some reasonable 
approximations (DVCS leading-twist and leading-order dominance, neglect of some ^ 
terms in the analytical expressions, ...), these set of relations can be linear. Eqs. [27] 
give four examples of such relations. All others can be found in Ref. [55]. Then, if a 
quasi-complete set of DVCS experimental observables can be measured at a given {xb, 
Q^, —t) point, one can build a system of eight linear equations with eight unknowns, 
i.e. the eight CFFs. 

Such system can be solved rather straightforwardly with standard matrix inversion 
and covariance error propagation techniques. This approach has been applied in 
Ref. |112] to the HERMES data. We recall that HERMES has this unique characteristic 
to have measured all beam-target single- and double-spin DVCS observables. The 
absence of cross section measurement at HERMES means however that these spin 
observables are actually under the form of asymmetries, i.e. a ratio of polarized cross 
sections to unpolarized cross sections of the form This has the consequence that 
the mapping is not fully linear and that some additionnal (reasonable) approximations 
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Figure 39. The Him CFF at —t ~ 0.28 GeV^ as a function oixB- Tlie empty squares 
sliow tlie result of tlie 7 CFFs free parameters fit from Refs. [IMl |T27l [1281 [T29]. The 
point at xb ~ 0.09 is from the fit of the HERMES data, the ones dX Xb ~ 0.25 and 0.35 
from the fit of the CLAS data and the one aX Xb ~ 0.36 from the fit of the JLab Hall 
A data. The solid red curves show the result of the VGG model (with bvai = &sea=l 
and without any D-term for H) . The solid dashed blue curves show the results of the 
model-based fit of Ref. [113) . The dash-dotted green curve shows the result of Him at 
-t = Q GeV^, i.e. the PDF (MRST02). 



have to be made such as the dominance of the BH squared amphtude over the DVCS 
squared amphtude. 

With an appropriate selection (with some partial redefinition) of eight HERMES 
observables, the others serving as consistency checks, the authors of Ref. |112] has 
been able to solve the system of eight equations and extract the eight CFFs with their 
uncertaintities. In this well-constrained approach, the uncertaintities on the CFFs refiect 
essentially the errors of the experimental data. Fig. 138) shows the results of this mapping 
technique for the three CFFs Him-, Him and H^e in the HERMES column (black circles). 
With the precision of the HERMES data, only the Him CFF come out to be clearly 
different from zero, all others CFFs being compatible with zero within error bars. 

In Fig. |38l the agreement of the mapping technique with the "brute force" least- 



GPDs in the valence region from DVCS 



81 



square minimization technique discussed in the previous section is striking (we note that 
the least-square minimization technique was also applied in Ref . |112] , with results well 
in agreement with those of Ref. |128] ). We refer the reader to Ref. [112] for a discussion 
on the reason why the error bar on Him for the third t value (around -.2 GeV^) is 
somewhat larger than for the other t values. We also note that HERMES measured 
the AuL asymmetry at four t values which should therefore allow in principle to extract 
Him at these four t values. However, the extracted Him CFF at the lowest | t \ point 
(around -.03 GeV^) turns out to be negative (although compatible with zero within 
two standard deviations, which does not discard a statistical fluctuation effect). This is 
both in the present mapping approach and in the least-square minimization approach 
discussed in the previous section. This negative value for Him at very low | t | is a bit 
surprising as, unless skewdness effects introduce a sign flip, it would imply a negative 
proton polarized parton distribution function. We do not display it here but it is shown 
in Ref. [TT2] . 

4-2.3. Fitting with only H One limitation of the two previous methods that we just 
presented is that every {Ee,XB,Q'^,t) kinematic point is taken individually and fitted 
independently of all others, in particular of its neighbors which have no influence on 
the fit. On the considered experimental data sets, it turns out that the resulting 
CFFs display a rather smooth behavior as a function of —t and don't show oscillations. 
This, in a way, vahdates the method but nothing prevents the occurrence of oscillations 
when studying other measurement^. Still, one could wish to improve the procedure 
and, staying in an almost model-independent fitting framework, enforce some general 
properties like the smoothness or continuity of the CFFs or the implementation of 
dispersion relations. 

One attempt to enforce the smoothness of the CFFs was made in 2009 in Ref. |125j . 
working with the CLAS beam spin asymmetries and the JLab Hall A unpolarized and 
beam-polarized cross sections and assuming the dominance of the GPD H. As we saw 
in section 14. 1[ this assumption is supported by all models and is expected to work at 
the 20 to 50 % accuracy. In Ref. [125] , this influence of the other GPDs was probed by 
fitting the data assuming that E, H and E either vanish or take their VGG values, and 
studying the dispersion of the fit results. 

To enforce the smoothness, one imposes a generic functional form to the H GPD. 
In Ref. [125j . the singlet combination H^ is parameterized in the dual model framework 
according to Eq. fl83l) (the GPD E is neglected). The t-dependence of the Bni coefficients 
is parameterized as: 

with to = - 0.28 GeV^. Such a parameterization correlates the x and t dependences. 
The reference scale is defined by Qo = 3 GeV^ and the evolution is performed with 

* Such an oscillating behavior is expected when studying only a single type of measurements, since the 
previous CFF fitting approach leads to an underconstrained problem. 
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Figure 40. Q^-behavior (1 < < 4 GeV^) of the extracted values of Hirn and 
of local fits (left) and global fit (right) on Hah B kinematics : 0.09 < -t < 0.2 GeV^ 
(a), 0.2 <-t< 0.4 GeV^ (b), 0.4 < ~t < 0.6 GeV^ (c), and 0.6 < -t < 1. GeV^ (d). 
The error bars include both statistics and systematics. iJ/„i ranges between and 10, 
Hfie between -7.5 and +7.5. Note the change of notational conventions with respect to 
Eq. p5|) and Eq. ^T9\\ to match the results published in ref. (I25j . The black full circles 
correspond to a;_B=0.125, red squares to 2:5=0. 175, green up triangles to 2:3=0. 250, 
blue down triangles to xb=0.360 and magenta open circles to 2:^=0.491. 



3 active quark flavors and A = 373 MeV. 

The convergent series Eq. fl87|) is truncated at some maximal value of A^max- A 
rough estimate of the uncertainty related to the specific choice of the truncation was 
obtained by comparing fit results with different values of A'max- Choosing Amax = 2, 
3 or 4 is sufficient to obtain reasonable fits to the data. Larger values of A'max produce 
numerical instability with some coefficients left largely undetermined; the overall quality 
of the fit also becomes poor. 

The fits to the Hall A and CLAS data were performed in two ways: either "locally" , 
i.e. fitting each individual kinematic bin {Ee,XB,Q'^,t) independently of the other, as 
described in previous section, or "globally", i.e. fitting all {Ee,XB,Q'^,t) kinematic bins 
simultaneously. The left and right panels of Fig. HO] display the results for the Hjm and 
Hjif. CFFs respectively, for both methods. The results for both kinds of fits are almost 
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always compatible, which is a good consistency check. As expected, the results of the 
global fits are in general smoother, due to the implementaton of the functional form 
for H. This is especially true concerning H^e '■ in the case of the local fits, this CFF 
shows large fluctuations between neighboring bins, with in some cases values falling even 
outside the plot range. The results of the local fits suffer from large fluctuations as the 
fits are not constrained enough in some bins but have the advantage of being almost 
mo del-indep endent . 

Both local and global fits give results with comparable accuracy for Hjm- While 
i^/m is rather precisely extracted, if^e is still poorly known, which was also a conclusion 
of section 14.2.11 

4 -2. 4- Neural networks For the sake of completeness, we finally mention the pioneering 
work Ref. |130j which constitutes the first attempt to extract CFFs from DVCS data 
with neural networks instead of traditional least square minimization methods. We 
will not cover it in details since it deals only with a subset of existing data, namely 
HERMES beam spin and beam charge asymmetries. This approach yields uncertainties 
in agreement with those obtained from fits with an a priori given functional form and 
standard statistical procedures. Being largely model- independent, the uncertainties 
estimated when extrapolating to the t — )• region are presumably safer. This body of 
results will certainly trigger new studies in the future. 

5. The future 

So far, the data relevant to the DVCS physics in the valence region have come from 
the JLab Hall A and CLAS experiments using a 6 GeV electron beam and from the 
HERMES experiment using a 27 GeV electron or positron beam. As described in detail 
above. Hall A has measured unpolarized and beam-polarized cross sections [58], CLAS 
beam spin asymmetries and longitudinally target spin asymmetries [61] and HERMES 
the complete set of beam charge, beam spin and target spin asymmetries [611 ESI [661 EZl 
EH EHl EQI |71]. We have made use of practically all these existing data in the previous 
sections by comparing them to model calculations or fitting them in order to extract 
CFFs. These data have allowed to show the successes and the limits of the present GPD 
parameterizations and to develop the first GPD of CFF fitter codes. Putting together 
all these informations, one can consider that the Hjm CFF is relatively well constrained 
and known at the ~15% level and that some first constraints on the Him and H^e CFFs 
start to appear. 

Although possibly more information can still be extracted from these data, one 
clearly wishes to have, ideally, more observables, more precise data and a larger phase 
space coverage. In the intermediate xb region 0.1), unfortunately not much more 
can be expected from HERMES since the experiment has shut down a few years ago. 
However, the COMPASS experiment with a 160 GeV muon beam is scheduled to have 
around 2016 a dedicated DVCS program with a specific recoil detector to ensure the 
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exclusivity of the process |131] . It should then be able to explore the xb range between 
0.01 and 0.1, thus with some partial overlap with HERMES. 

On a shorter time scale, a lot of new data are expected to come from JLab. The 6 
GeV era has just finished in summer 2012 and many data are currently under analysis: 

• In Hall A, the experiment E07-007 [132] has carried out DVCS measurements 
at two beam energies (6 and 4 GeV) for which it is planned to extract the 
unpolarized and beam-polarized cross sections. At fixed xb and Q^, the beam 
energy dependence, analog to a Rosenbluth separation, will allow to separate the 
pure DVCS contribution from the BH-DVCS interference contribution. Strong 
contraints on the real CFFs, in particular -ff/je, are expected from this measurement. 

• With CLAS, the experiment E06-003 [133j has carried out DVCS measurements at 
~ 6 GeV with a polarized beam. A first set of beam spin asymmetries released 
of this experiment has already been published [61] and discussed in the previous 
sections but another set with double statistics is currently under analysis. Most 
importantly, unpolarized and beam-polarized cross sections are in the process of 
being extracted. A glimpse on these cross sections is available in Ref [60]. These 
cross sections are expected to be less precise than the Hall A ones due to the 
systematic uncertainties inherent to a large acceptance detector such as CLAS. 
However, they will cover a much larger phase space. Strong constraints on the H^e 
and Hirn are expected from this measurement. 

• With CLAS, the experiment EOS- 114 |134] has carried out DVCS measurements 
at ^ 6 GeV with a longitudinally polarized proton (and neutron) target and a 
polarized beam. Improved Aul (w.r.t. Ref. [63]) and, for the first time. All 
measurements are thus expected soon. This will allow to further constrain the Him 
CFF in particular. 

JLab is currently in an upgrade phase and plans to deliver a 12 GeV beam around 
2015. A DVCS program in Hall A and in CLAS has already been approved: 

• In Hall A, the experiment E12-06-114 |135] will measure the unpolarized and beam- 
polarized cross sections in a new kinematical regime (smaller xb and larger Q"^). 

• With CLAS, the experiment E12-06-119 |136] will use a polarized beam and a 
longitudinally polarized target to measure unpolarized cross sections and beam 
spin, target spin and double spin asymmetries. 

• With CLAS, the proposal C12-12-010 |137] aims at measuring the DVCS reaction 
with a tranversally polarized target. This should bring strong constraints on the E 
GPD (see Eq. (jSO])). 

• Besides the DVCS experiments off the proton, the experiment E12- 12-003 |138] 
aims at measuring DVCS on the neutron. Except for the pioneering measurement 
of the Hall A [59], DVCS on the neutron has never been measured. It is obviously 
an indispensable process to measure in order to perform a flavor separation of the 
GPDs. 
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Fig. HI] compares the {xbiQ"^) domains that are explored or will be explored by 
JLab 12 GeV, HERMES, COMPASS and Hl/ZEUS regarding the DVCS and DVMP 
processes. This illustrates the complementarity of all these facilities, the near future 
belonging to the JLab 12 GeV and COMPASS facilities. In the following two sections, 
we show some examples of what is expected to be achieved with the JLab 12 GeV facility 
and in the third one a comparison of the predictions for the various models presented 
and discussed in sections |3] and H] for COMPASS kinematics. 
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Figure 41. (xb.Q'^) domain explored by JLab 12 GeV, HERMES and COMPASS 
regarding tlie DVCS (the limits take into a. W >2 GeV cut and an approximate 
estimation of the luminosity of each facility). 

5.1. Hall A 

While the CLAS12 detector will explore a wide phase space region for the DVCS process, 
the DVCS program in Hall A will be to focus on some specific kinematics and make 
precision measurements. In terms of systematic uncertainties, we recall that the typical 
momentum resolution of the Hall A arm spectometers is of the order of 10~^ (compared 
to 10^^ for CLAS12) and that, in terms of statistics, the luminosity that can be reached 
in Hall A is of the order of 10^^cm~^s~^ (compared to 10^^cm~^s~^ for CLAS12). In 
particular, before thinking of extracting GPDs or CFFs out of DVCS data, it is of the 
utmost importance to ensure that the "handbag" formalism is applicable, in particular 
at the relatively low values that can be reached at JLab. One of the signatures to be 
looked for is the scaling behavior of the CFFs, i. e. the property that they don't depend 
on at leading-twist. 

The preliminary tests of scaling carried out at 6 GeV by the Hall A collaboration 
are encouraging (see Fig. [T3l) but are limited in the range (between 1.4 and 2.4 
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GeV^). Fig. WI\ shows the gain in the range that can be obtained with the JLab 
12 GeV beam energy increase. The errors on this figure were estimated for 90 days of 
beam time and the running conditions of the JLab Hall A proposal [135] . 
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Figure 42. Projected uncertainties from the JLab 12 GeV Hall A E12-06-114 
experiment [135) for the Q^-dependence of the Him CFF for different x b (figure done 
by C. Munoz). 



5.2. CLAS12 

CLAS12 is expected to measure all the DVCS observables accessible with a polarized 
beam, a longitudinally and a transversely polarized target. Except for beam charge 
observables, this makes up for a complete program. In this section, we show what 
can be achieved in terms of the extraction of CFFs using the technique presented in 
section 14.2.11 This study has been done in collaboration with H. Avakian. 

For each (xs, Q^, t) bin, the distributions of the various independent DVCS 
spin observables which are measurable with a longitudinally polarized beam and a 
longitudinally or/and transversely polarized proton target, have been generated: A^u, 
^UL, All, Ajjx, Auy, A^^, A^y (in addition to the unpolarized cross section). The 
VGG values [H2] for the four GPDs H, E, H and E have been used to generate these 
distributions. Then, ep — )■ ep7 events generated according to the BH+DVCS cross 
sections have been fed into a fast Monte-Carlo simulating the CLAS12 acceptance and 
efficiency. Assuming 80 days of beam time for the unpolarized target run at a luminosity 
of 10^^ mc~^s~^ (from which the unpolarized cross section and Alu are planned to 
be extracted), 100 days of beam time for the longitudinally polarized target run at a 
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Figure 43. Resulting i/^e CFF from the simultaneous fit of A^c/, Aul, All, Au^: 
Auy, Alx, ALy and of the unpolarized cross section, for each {xb, Q^, t) bin with the 
fitting code of Refs. [1261 HSU II2S 1129]. The extracted CFFs for which the error bar 
was larger than 3 were removed, (study done in collaboration with H. Avakian) . 

luminosity of 2.10^^ mc^^s^^ (from which and An are planned to be extracted), 
100 days of beam time for the transversely polarized target run at a luminosity of 
5.10^^ mc^^s^^ (from which Au^, Auy, A^^ and A^y are planned to be extracted) and 
furthermore assuming 80% target polarization, one can assign (statistical) error bars to 
the distributions and then fit them. The goal is to extract the seven CFFs : if_Re, 
Ejie, Hjie, Ejie, Hj^n, Ej^n, Hj^. As already mentioned in section H. 2.11 Ejm=0 is set to 
in this study. 

The results for the seven CFFs issued from the simultaneous fitting of the 
distribution of A^u, Am, ^ll, ^c/x, ^r/y, ^Lx, A^y and of the unpolarized cross section 
with the procedure of Refs. [1261 [EH] are displayed, for all {xb, Q^ t) bins 

in Figs. Us] to HHl The reconstructed CFFs with error bars should be compared to the 
generated ones which are represented by the solid curve. The panels where there are 
blue solid curves and not CFF reconstructed means that the particular fitting code that 



GPDs in the valence region from DVCS 



88 




Figure 44. Resulting Him CFF from tlie simultaneous fit of Alu, Aul, All, Au^, 
Auy, Alx, ALy and of the unpolarized cross section, for each {xb, Q^, t) bin with the 
fitting code of Refs. [1261 HSU II2S 1129]. The extracted CFFs for which the error bar 
was larger than 150% were removed, (study done in collaboration with H. Avakian). 

we have used was not able to reconstruct the CFFs rehably. In Figs. US] to SHI we have 
indeed not plotted the extracted CFFs for which the error bar was too large (see the 
captions of the figures was the detailed criteria). This does not mean that there are no 
data reconstructed in those bins. Other fitting algorithms (probably with some model- 
dependent inputs) can certainly make use of these data and extract some constraints 
on the CFFs and GPDs. 

The dashed curves on the imaginary parts of the CFF plots show the corresponding 
values of the GPDs for zero skewness argument, i.e. H{x, 0, t), H{x, 0, t) and E{x, 0, t), 
according to VGG. These latter are the quantities which have a simple probability 
interpretation. One therefore sees from these (model-dependent) curves, showing the 
difference between H{x,x,t) and H{x,0,t), the effect of the skewness. We will come 
back to this issue in the section [6l 

With such complete experiments, which comprise all DVCS observables except 
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Figure 45. Resulting Hni. CFF from the simultaneous fit of A^u, Ajjl, All, Ajj^, 
Auy, Alx, A^y and of the unpolarized cross section, for each {xb, Q^, t) bin with the 
fitting code of Refs. [1261 HSU 1121 |l29]. The extracted CFFs for which the error bar 
was larger than 0.3 were removed, (study done in collaboration with H. Avakian). 

for beam charge asymmetry, we observe that the seven CFFs can be reconstructed for 
essentially all {xb, Q^, t) bins with quite good precision. In particular, the measurement 
of the transverse target spin observables is crucial to reconstruct the CFFs related to the 
GPD E. The same study without such observables allows the reconstruction of the CFFs 
related to the GPDs H and H only, leaving the CFF Ejm practically unconstrained. 

In Sec. 12.2! and Sec. 14.1! we already insisted on the required accuracy in the 
determination of the kinematics of observables. The BH cross section varies strongly 
when the outgoing photon is emitted in the direction of the incoming or scattering 
electron. An variation of 1 % in a;^ can induce a variation of 10 % in the cross section. 
This has some important consequences regarding the CFF fitting method described 
above. For the CFFs to be accurately reconstructed, it is necessary to have a complete 
set of observables on the same {xb, Q^, t) bin. If this is not the case, we have to work 
with neighboring bins (for example with similar, but not equal xb)', this approximation 




Figure 46. Resulting Hjm CFF from tlie simultaneous fit of Alu, Aul, All, Au^, 
Auy, Alx, ALy and of the unpolarized cross section, for each {xb, Q^, t) bin with the 
fitting code of Refs. [1261 HSU 1121 |l29]. The extracted CFFs for which the error bar 
was larger than 150% were removed, (study done in collaboration with H. Avakian). 

will generate a systematic uncertainty on the reconstruction of CFFs that can be quite 
large. 

5.3. COMPASS 

The COMPASS experiment plans to measure the correlated beam charge-spin 
observables: 

Scs,u = (T^^ + (^'^ , (135) 

Dcs,u = <y^^ - <y-^ 

since the muons of the beam originate from pions decays which induces a correlation 
between the charge and the spin of the muons. 
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Figure 47. Resulting iS^je CFF from tlie simultaneous fit of Aijl, All, Ajj^, 

Auy, Alx, A^y and of the unpolarized cross section, for each {xb, Q^, t) bin with the 
fitting code of Refs. [1261 HSU 1121 |l29]. The extracted CFFs for which the error bar 
was larger than 1.5 were removed, (study done in collaboration with H. Avakian). 



In Fig. inni we show the predictions of the four models which we discussed 
in sections [3] and HI i.e. the VGG, GK, dual parameterization and KM models. 
Calculations have been done for ii^^=160 GeV, xb=0.05, Q^=2 GeV^ and —t=-0.2 
GeV^. All four models show relatively similar features. The differences lie in the global 
normalization of the Scs,u observable and in the behavior at the lowest and largest 
values of the Dcs,u observable. 
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Figure 48. Resulting Eim CFF from tlie simultaneous fit of Ajjl, All, Ajj^, 

Auy, Alx, A^y and of the unpolarized cross section, for each {xb, Q^, t) bin with the 
fitting code of Refs. [1261 HSU 1121 |l29]. The extracted CFFs for which the error bar 
was larger than 150% were removed, (study done in collaboration with H. Avakian). 
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Figure 49. Resulting E^e CFF from the simultaneous fit of Am, Aul, All, Ajj^, 
Auy, Alx, A^y and of the unpolarized cross section, for each {xb, Q^, t) bin with 
the fitting code of Refs. [1261 HSZl USS [US], (study done in collaboration with H. 
Avakian) . 
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Figure 50. Predictions at E^,^im GeV, xb=0.05, Q'^=2 GeV^ and -^=-0.2 GeV^ for 
tire 0-dependence of Scs,u (left panel) and Dcs,u (right panel) as defined in Eq. 11351 
Solid red curve: VGG; dashed blue line: GK; dotted black: dual model; dot-dashed 
green: KM model. 
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6. Prom Compton Form Factors to spatial densities 

As described above, the GPDs at ,^ = 0, e.g. H{x,0,t), are mapping out the combined 
probabilities in transverse position and longitudinal momentum of the quarks in the 
nucleon (see Eq. ([8])). We saw in the previous section that CLAS12 will allow to 
essentially extract all CFFs, with more or less precision depending on the OFF and 
the kinematics, over the range 0.1 < < 0.6 and t^m ^ ~^ ^ 1 GeV^. In particular, 
if we focus on the unpolarized GPD H, the CFF Hjm = H{^, ^, t) — H{—^, ^, t) can be 
extracted quite precisely. In the following, in this pioneering exercise, we will make the 
approximation of neglecting the antiquark contribution to the Him CFF, i.e. neglect 
H{—^,^,t) w.r.t. H{^,^,t). At CLAS kinematics, according to the GK and VGG 
models, H{—^,C,,t) is about 20% of Hjm while at HERMES kinematics, it is about 
30%. This approximation being clearly set, given the uncertainties on H{C,,C,,t), and 
modulo a (model-dependent) skewness correction of the form H{C^,0,t)/H{C,,C,,t), one 
can adress several questions : 

• With which accuracy can one extract H{x, b±) from the measurement of the 



• How does one perform such error propagation ? 

• What is the model dependence of the skewness correction ? 

Eq. ([8]) can be equivalently expressed, for a circularly symmetric function, through 
the Hankel transform: 



where A_|_ = |A_|_|, and Jq is the Bessel function of order 0. 

We present here a simple numerical algorithm which adresses the error propagation 
in this transform. We illustrate it by focusing on Him and by taking one particular 
{xByQ"^) bin in Fig. HU e.g. (0.1,2.5 GeV^). Fig. [5l] (left panel) displays the pseudo- 
data and the associated errors contained in this bin. The procedure consists in smearing 
"vertically" the seven values of H{C,, ^, t) (which correspond to seven t values) according 
to a Gaussian distribution with a standard deviation equal to the error bar of the point. 
These 7 "new" points are then fitted by a function, which we take as an exponential 
Ae^^ with two free parameters, the normalization A and the slope B. At fixed a;^, an 
exponential ansatz is motivated by most GPD models (for instance VGG, GK, KM,... 
which we discussed in section [3]). In principle, any other fit function can be used, and 
can serve as a way to estimate a systematic uncertainty associated to this method. This 
procedure (smearing + fitting) is repeated several thousand times so that one obtains 
several thousand exponential functions, shown in Fig. [5T] (middle panel). 

Then, each of these several thousand exponentials is transformed through Eq. (11361) 
so that one obtains a series of Hankel transforms, now as a function of bj_. This transform 
can be done analytically in the present case but any numerical method is also possible 
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Figure 51. Left: pseudo-data for the t-distributfon of the Him CFF (zoom of the 
{xb,Q^) bm (0.1,2.5 GeV^) of Fig. HH). Middle: same as left figure with, superposed, 
the series of Ae^* red curves fitting the smeared pseudo-data. The solid black curve 
shows the VGG Him CFF, i.e. H{£_,^,t). The dashed black curve shows the VGG 
value for H{£^,Q,t). Right: same as left figure with, superposed, the series of Ae^* 
red curves fitting the smeared pseudo-data, where these latter have been corrected by 
HitO,t)/Hi^,^,t). 

if the function fitting if ^, t) does not have a simple Hankel transform. The idea is 
then to look, for various values of b±, at the dispersion of all the Hankel transforms. 
The spread of the Hankel transforms for a few (arbitrarily) selected b± values is shown 
in Fig. [521 One sees that the resulting distributions are quasi-Gaussian. One can 
therefore fit those distributions by a Gaussian function from which one extracts the 
mean and the standard deviation. The black points in Fig. [53] show the result of this 
procedure. The error bars at each b± value result from the propagation of the error 
bars of H{C,,C,,t) displayed in Fig. [SH We choose to display here only seven b± points 
but the procedure can be applied to any number of b± points so that one can obtain a 
continuous distribution as a function of b±. The interest of this method is to properly 
propagate the errors on H{C,, ^, t) to its Hankel transform. In particular, it takes into 
account the correlations between the parameters used to fit H{^, ^, t) (in the present 
case, the normalization and the slope of the exponential). 

However, in order to be able to interpret such distribution as the transverse spatial 
density of the quarks at the particular value of x = ^ = 0.053 taken here, a correction 
needs to be made. So far, we have Hankel-transformed H{^,^,t). However, the spatial 
density interpretation requires the knowledge of H{^,0,t). The correction to pass 
from one to the other has to be model-dependent since the x and ^-dependencies 
cannot be measured independently with the DVCS process. A model-dependent 
H{^,0,t)/H{^,^,t) skewness correction factor should therefore be applied. 

Fig. [5l] shows such deskewing correction factor for three models that we discussed 
in section [3l It is seen that below —t=l GeV^, they don't differ by more than 10%. 
They all three correspond to a correction of the order of 20% at small t with a similar 




GPDs in the valence region from DVCS 



97 



100 

80 
60 
40 
20 


100 
80 
60 
40 
20 


100 
80 
60 
40 
20 






1.5 



2 0.5 



5 2.5 




0.2 



0.4 







_ b^=1 .75 / L\ 

n \ 

' \ 




/ 1 

V \ 








1 " ...Tfl _ , _ , 


; / \ 

/ \\ 
'Jl , , , > 


1 „ , 


0.1 


0.02 


0.04 



Figure 52. Distributions of the Hankel transforms H{x — 0.053, for selected 6j_ 
values, for the series of fit curves corresponding to the middle panel of Fig. [51] The 
mean and standard deviation (respectively the parameters P2 and P3 in the inserts 
of each panel) of the Gaussian fit of these distributions is extracted to obtain the 
&_L-dependence of the spatial density. 



evolution as —t grows up to 1 GeV^. Such comparison between different models can 
serve to estimate a systematic uncertainty in the deskewing. We note that the deskewing 
factor depends on t in a similar way for all three models. As | t | grows, the deskewing 
factor grows towards one (and can even go over one). This means that it flattens the t 
slope of H{C,,0,t) w.r.t. to the measured or extracted H{^,^,t). All three models also 
show a similar x^-dependence between the CLAS and the HERMES kinematics, i.e. 
the deskewing correction is less important as increases. 

We then apply such a deskewing factor to the pseudo-data of our selected {xb, Q^) 
bin. This yields the new set of data shown in Fig. |5l] (right panel). Then, the procedure 
that we just described is applied to these "deskewed" data (see the red curves in Fig.lSTj 
right panel). This results in the the h±_ distribution given by the red points in Fig. [53| 
which can now be properly interpreted as a spatial density. We stress that we place 
ourselves in a leading twist and leading order framework in this whole exercise. 
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Figure 53. 6^-dependence of the spatial density H{x = 0.053, 6^) resulting from the 
Hankel transform of Fig. 1511 middle panel, i.e. without any "deskewing" correction 
(black points), and of Fig. [511 right panel, i.e. with "deskewing" correction (red 
points). 



We have focused so far on one particular [xbiQ"^) bin in order to illustrate the 
method. We now process a series of bins in order to provide an imaging of the 
nucleon. We select in Fig. HH the Him row at Q^=2.'b GeV^ where there are 7 xb 
bins (if corrections are under control or found negligible, different rows can of 
course be combined). Fig. 155) shows the result where we now see the evolution of the 
transverse spatial density as a function of x^. In particular, one notes that, as 
decreases the radius of the proton increases. This is a feature which was implemented 
in VGG and it serves as a proof of principle that, after the different steps, of computing 
the DVCS observables from the VGG model, generating the pseudo-data according 
to these, processing these through a simulation of the CLAS12 detector, fitting these 
pseudo-data in order to extract the Him CFF and the final Hankel-transform procedure 
that we just discussed , one recovers the original features of the VGG model. 

We finally apply this procedure to real data. We recall that Him has been extracted 
out of the CLAS and HERMES data (see Fig. [38]) . Figs. [56] and 1571 show the results of 
the procedure. Of course, the Him data have much larger error bars than the CLAS12 
simulations and the exponential fits (red curves on the left panels of Figs. [56] and [571) 
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Figure 54. "Deskewing" factor H{^,0,t)/H{^,^,t) as a function of -t for the VGG 
model (red curves), the GK model (blue curves) and the dual model (black curves). 
The solid curves correspond to a;B=0.1 (HERMES kinematics) and the dashed ones to 
xb=0.25 (CLAS kinematics). 

have a lot of dispersion. The distribution of the Hankel transforms at fixed b± are then 
not as clean as in Fig. [521 In particular, the fits in Figs. |56] and |57] sometimes are 
horizontal straight lines (which correspond to exponentials with a zero slope) and this 
produces double peaks in the b± slices. Nevertheless, a dominant Gaussian-like peak 
always remains distinguishable from which a centroid and a standard deviation can be 
extracted, yielding the spatial charge densities of Fig. [56] (right panel), corresponding 
to xb=0.25, and Fig. [57] (right panel), corresponding to 0:^=0.09. Given the size of the 
experimental errors on the Hj^ data, we choose to ignore in this exercise the uncertainty 
associated to the deskewing. We also neglect at this stage the uncertainty associated 
with the fact that the different DVCS-BH asymmetries that are used to extract Hj^ are 
given at similar, but not strictly equal, kinematic {xBitiQ"^) points. This being said, in 
spite of the large error bars and some uncertainties in the procedure that we ignored, 
we can distinguish the same features than we observed with the simulations, namely, 
an increase of the size of the proton with xb decreasing. The overall increase in the 
normalization of the density is due to the rise of the unpolarized parton distribution 
function as xb decreases. As a further illustration, we fitted the spatial charge densities 
of the right panels of Figs. [56] and [57] by a Gaussian function which we plot as contour 
plots in Fig. [58] 
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Figure 55. Left: Evolution of ttie &j_-spatial densities as a function oi xb- The seven 
curves correspond to the Hankel transforms of the Q^=2.b GeV^ row of Him in Fig- HH 
Right; alternative view in color for visibility. 
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Figure 56. Left: the black squares show iJ/„i extracted without any deskewing 
correction and the black stars with deskewing correction, from the CLAS data. For 
visibility, the t values of the squares and stars are slightly shifted around the true t 
values. The red curves are the fitted Ae^^ functions obtained from the smearing of 
the Him points. Right: resulting spatial charge density. The blue dashed curve is the 
result of a fit by a Gaussian function. 
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Figure 57. Left: the black squares show Hi„i extracted without any deskewing 
correction and the black stars with deskewing correction, from the HERMES data. 
For visibility, the t values of the squares and stars are slightly shifted around the true 
t values. The red curves are the fitted Ae^^ functions obtained from the smearing of 
the Hi„i points. Right: resulting spatial charge density. The blue dashed curve is the 
result of a fit by a Gaussian function. 




Figure 58. Left: contour plot of the spatial charge density at CLAS kinematics 
(a;_B=0.25) extracted from Fig. [56l Right: contour plot of the spatial charge density 
at HERMES kinematics (a;B=0.09) extracted from Fig.[57l 
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7. Conclusions and outlook 

In this work, we have briefly reviewed the fleld of the Generalized Parton Distributions 
and Deeply Virtual Compton Scattering in the valence region, with emphasis on 
the information which can be extracted from present and forthcoming data. After 
briefly recalling the theoretical formalism and the properties and interests of the 
GPDs, we have reviewed the few existing data from the HERMES and JLab 
facilities. Then, wc presented four widely used GPD parameterizations based on 
different approaches (double distributions, dual parameterization, and Mellin-Barnes 
representation), followed by a more general discussion of DVCS observables within a 
dispersion relation framework. In the present work, our study was based on a leading- 
twist and leading-order QCD assumption. We compared the results of these different 
approaches to the existing data. We found that, although many features and trends 
of the experimental observables were correctly accounted for, no model satisfactorily 
describes simultaneously the full set of current HERMES and JLab data. We also 
presented an alternative approach of fitting the Compton Form Factors, which are 
functions of GPDs, in more or less model-independent ways. Given the few existing 
data, the information that can be extracted in this way is limited but nevertheless 
brings additional clues on the behavior of GPDs. 

Putting all the model-dependent and model-independent pieces of the puzzle 
together, one can conclude so far that the CFF Hi^ can be considered as known in 
the valence region at the ^ 15% level and, to a lesser extent, the CFF Him- We have 
further shown that an imaging of the nucleon starts to appear from these analyses. 
One feature which emerges from this analysis is that the increase of the t-slope of 
Hf„i with decreasing is reflecting the increasing transverse size of the nucleon as one 
probes partons with smaller and smaller momentum fractions (in the so-called light-front 
frame). This yields an image of the nucleon with a core of valence quarks surrounded 
by a cloud of quark- antiquarks. Another feature, resulting from the observation that 
the CFF Him has a flatter ^-dependence than the Him and barely varies with is that 
the axial charge of the nucleon tends to stay concentrated in the core of the nucleon. 

Furthermore, we explored the future of the fleld, in particular by summarizing all 
the simulation work that has been carried out for the JLab 12 GeV facility. With 
the upgrade in energy and luminosity (for CLAS) of JLab, an unprecedented set of 
data in precision and phase-space coverage is anticipated. We showed that essentially 
all leading-twist CFFs will be extracted, some with better precision than others. In 
the flnal chapter of this review, we flnally proposed a method to transform a CFF 
measurement into a spatial charge density measurement with proper error propagation. 
We apphed this technique to the JLab 12 GeV pseudo-data as well as to the few existing 
data from JLab 6 GeV and HERMES. 

The future of the fleld is bright with numerous new data in perspective from JLab 
12 GeV and COMPASS. The few existing data, with their limitations in terms of 
precision and phase space coverage, allow to develop the techniques (models, flts,...) to 
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be employed in the future. Within the present work, we showed within the leading-twist 
and leading order QCD assumptions, a proof of principle to extract the GPD related 
observables from the deeply virtual Compton process. Furthermore, on can anticipate 
further developments on the theory side in calculating the corrections required by the 
anticipated precision of the forthcoming data. Such program will allow to image the 
partonic structure of the nucleon. Visualizing will then deepen our understanding. 
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